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Abstract 

We present the Monte Carlo event generator YFSWW3 version 1 . 16 for the pro- 



Q I cess of I^-pair production and decay in electron-positron collisions. It includes 

' C(ct) electroweak radiative corrections in the WW production stage together with 

the 0{a^) initial-state-radiation (ISR) corrections in the leading-logarithmic (LL) 
^ I approximation, implemented within the Yennie-Frautschi-Suura (YFS) exclusive 

^ ' exponentiation framework. The photon radiation in the W decays is generated by 

^ the dedicated program PHOTOS up to 0{a^) LL, normalized to the W branching 

ratios. The program is interfaced with the r decay library TAUOLA and the quark 
fragmentation/hadronization package JETSET. The semi-analytical code KorWan for 
the calculations of the differential and total cross-sections at the Born level and in 
the ISR approximation is included. 

To be submitted to Computer Physics Communications 



f Work partly supported by the Polish Government grant KBN 5P03B12420, the European Commis- 
sion 5-th framework contract HPRN-CT-2000-00149, the US DoE Contracts DE-FG05-91ER40627 
and DE-AC03-76ER00515, and the Polish-French Collaboration within IN2P3 through LAPP An- 
necy. 



CERN-TH /2001-017 
DESY 01-027 
January 2001 



Contents 

I Introduction 



2 Physics Contents of YFSWW3 

|2.1 Notationl 

2.2 Master .bbrmula 

12.2.1 YFS Infrared Factori 

2.2.2 Coulomb Correction 

2.2.3 Anomalous Triple- Gauge Couplings . . 

2.2.4 (j's for H/-Pair Production| 

2.2.5 Pretabulated EW Corrections 

2.2.6 Fixed-Order Calculation 

2.2.7 Radiative Corrections for W Decays] . . 

2.2.8 ISR Reference Differential Distribution 

2.3 YFS Reduction/Extrapolation Procedures . . 

2.4 Multichannel YFS Reduction/Extrapolation . 



3 Monte Carlo Algorithm 

3.1 Main Steps in the Construction of the MC Weight 

3.2 MC Weights and Absolute Normalization 



4 Structure of the Program 

4.1 demo - Demonstration Program . 

4.2 data_f iles - Default Input Data 

4.3 yf sww - Master Unit 



4.4 model - Matrix Elements 




4.5 ewe - Electroweak Library 



4.6 interf s - Interfaces to External Libraries . . 
O wdeclib - External Libraries for W Decay . . 

4.8 glib - Histogramming and Numerical Library 

4.9 semian - Semi-Analytical Program KorWan . . 

4.10 f f - FF Package! 

4.11 dok - Related Papers 

4.12 rewt - Reweighting 'Pools 

[4.12.1 Reweighting KoralW Events by YFSWW3| 
TTITl Reweighting YFSWW3 Events by YFSWW3| 
iJTS Reweighting YFSWW3 Events by KoralWl 



5 How to use the Program 






5.1 Principal Entries of YFSWW3| 


5.2 Input/Output Parameters 




5.3 Printouts ot the Program 



5 Summary 



A Program Parameters and Their Settings] 
B Output of the Demo Program 



35 
37 
44 



2 



PROGRAM SUMMARY 



Title of the program: YFSWW3, version 1.16. 

Computer: any with the FORTRAN 77 compiler and the UNIX/Linux operating system 
Operating system: UNIX/Linux 
Programming language used: FORTRAN 77 
High speed storage required: < 10 MB 

Size of distributed program: 3,237 kB 
Distribution format: tar gzip file (size of 1,103 kB) 

No. of cards in combined program and test deck: about 23,000 plus about 48,000 of 
auxiliary packages: KorWan, FF, JETSET, TAUOLA, PHOTOS. 

Keywords: 

Standard Model (SM), LEP2, linear colliders (LC), quantum electrodynamics (QED), 
quantum chromodynamics (QCD), boson W, VF-pair production, W decay, W branching 
ratio (BR), triple gauge boson couplings (TGC), quartic gauge boson couplings (QGC), 
four-fermion (4/) background, radiative corrections, Yennie-Frautschi-Suura (YFS) expo- 
nentiation, initial-state radiation (ISR), leading-log (LL) approximation. Coulomb effect, 
final-state radiation (FSR), electroweak (EW) corrections, leading-pole approximation 
(LPA), Monte Carlo (MC) simulation/generation. 

Nature of the physical problem: 

The process of the Vl^-pair production is important for precise tests of the Standard Model 
as well as searches for "new physics" at LEP2 and future linear coUiders. In order to match 
the experimental precision necessary for a successful physics programme, quantum effects 
(the so-called radiative corrections) have to be included into a theoretical description of 
this process. It turns out that not only the so-called universal corrections (initial-state 
radiation, the Coulomb effect, "naive" QCD corrections, etc.) are necessary, but also 
the 0{a) electroweak corrections in the WW production are needed to reach the desired 
theoretical accuracy. All these effects should, preferably, be included in a Monte Carlo 
event generator in order to account for realistic experimental set-ups. 

Method of solution: 

The Monte Carlo event generator for the combined W^-pair production and decay process 
including 0{a^) LL ISR effects, the Coulomb correction (usual or screened), the "naive" 
QCD effect, the 0{a) EW corrections in the WW production stage, implemented within 
the YFS exclusive exponentiation framework is provided. Multiphoton radiation in the 
WW production is generated according to the YFS MC method. The photon radiation in 
the W decays, normalized to the W BRs, is generated by the LL-type MC program PHOTOS 
(up to two photons). The decays of r's including radiative corrections are simulated 
by the dedicated package TAUOLA. The quark fragmentation/hadronization is performed 
with the help of the Lund program JETSET. The program can provide both weighted and 
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unweighted (weight = 1) events. Any experimental cut and apparatus efficiency may be 
introduced easily by rejecting some of the generated events. 

Restrictions on the complexity of the problem: 

The LPA is used (in two versions) to describe the signal WW production process. Mul- 

tiphoton radiation according to the YFS exponentiation scheme is generated only for the 
WW production stage. 0{a) EW corrections (in LPA) are included only in the WW 
production stage. The ISR effects beyond 0{a) arc included in the LL approximation. 
Non-factorizablc corrections (interferences between the production and decay stages) are 
approximated by the so-called screened Coulomb ansatz. Spin correlations between the 
WW production and decays are fully included only at the Born and the ISR levels. Ra- 
diative corrections in the W decays are included into an overall normalization through 
the W BRs, and the real photon radiation is generated in the LL approximation (up 
to two photons) by the program PHOTOS. Anomalous triple gauge boson couplings are 
included in the Born-like matrix element, i.e. with the universal SM corrections only. 
Quartic gauge boson couplings are implemented according to the Standard Model only 
(no anomalous couplings). The r decays and quark hadronization arc performed, respec- 
tively, with the help of the dedicated packages TAUOLA and JETSET. No 4/ background 
processes are included. 

Typical running time: 

200 CPU seconds of a PC Intel Pentium III @ 550MHz per 1000 unweighted events, for 
the parameter settings as they are given in the demonstration program. 
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1 Introduction 



The process of W^-pair production in electron-positron colliders is very important for 
testing the Standard Model (SM) and searching for signals of possible "new physics"; 
see e.g. Refs. One of the main goals when investigating this process at present and 

future e^e~ experiments is to measure precisely the basic properties of the W boson, such 
as its mass Mw and width Tw This process also allows for a study of the triple and 
quartic gauge boson couplings at the tree level, where small deviations from the subtle 
SM gauge cancellations can lead to significant effects on physical observables - these can 
be signals of "new physics". 

In this work we present a Monte Carlo (MC) event generator that simulates the produc- 
tion and decay of the W pair at colliders. The integrated cross sections and arbitrary 
differential distributions can be calculated from a series of constant-weight fully inclu- 
sive MC events. The program embodies the Standard Model scattering matrix element 
with the 0{a) radiative corrections for the doubly-resonant component of the scattering 
matrix element. The basics of the model used in the program and selected numerical 
results were already presented in earlier publications . The important components of 
the presented program/calculation are the complete 0{a) radiative correction for on-shell 
VF-pair production, which are implemented according toQ Refs. |]7|-p!0|, see also Ref. |TT| . 



The present YFSWW3 MC event generator was instrumental in achieving the new improved 
theoretical precision of 0.4% for the total cross section of the WW production process 
at the highest LEP2 energies, as a result of direct and detailed comparison with the 
RacoonWW MC program [13,0 within the 2000 LEP2 MC workshop (this conclu- 



sion was also supported by the comparison between RacoonWW and the calculation of 
Ref. HI), see Section 4.3 in Ref. ||. 

Let us briefly discuss the basic physics properties of the ly-pair production and de- 
cay process. Since the W^'s are unstable and short-lived particles, the W pairs are not 
observed directly in the experiments but through their decay products: four-fermion (4/) 
final states (which may then also decay, radiate gluons/photons, hadronize, etc.). As 
high energy charged particles are involved in the process, one can also observe energetic 
radiative photons. So, at the parton level, one has to consider a general process: 

e+ + e--^4/ + n7, (n = 0, 1, 2, . . . ), (1) 

where also some background (non-iyiy) processes contribute. In a theoretical description 
of this process - according to quantum field theory - one also has to include virtual effects, 
the so-called loop corrections. This general process is very complicated since it involves 
~ 80 different channels (4/ final states) with complex peaking behavior in multiparticle 
phase space and a large number of Feynman diagrams. Even in the massless-fermion ap- 
proximation, the number of Feynman graphs grows up from 9-56 per channel at the Born 



level to an enormous 3579-15948 at the one-loop level . The full one-loop calculations 

^ We are grateful to the authors of these works for providing us with the relevant parts of the computer 
code. 
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have not been finished yet, even for the simplest case (doubly plus singly H^- resonant dia- 
grams) |jl6[ . But even if they existed one would be faced with problems in their numerical 
evaluation in practical applications, particularly within Monte Carlo event generators - 
they would be prohibitively huge and slow. These are the reasons why efficient approxi- 
mations in the theoretical description of this process are necessary. These approximations 
should be such that on the one hand they would include all contributions/corrections that 
are necessary for the required theoretical accuracy (dependent on the experimental preci- 
sion) and on the other hand they would be efficient enough for numerical computations. 
Given the complicated topologies of the (4/ + ny) final states, such calculations should 
be, preferably, given in terms of a Monte Carlo event generator that would allow one to 



simulate the process directly [jT^|T^. 

Our solution to this consists of two complementary Monte Carlo event generators: 
YFSWW3 and KoralW. The latter includes the full lowest-order e~^e~ — > 4/ process, but 
with simplified radiative corrections - the universal ones such as initial-state radiation 
(ISR), the Coulomb effect, etc. In YFSWW3, on the other hand, the lowest-order process 
is simplified - only the doubly ly-resonant contributions are taken into account, but 
inclusion of the radiative corrections in this process goes beyond the universal ones. In the 
current version of YFSWW3 only those non-universal (non-leading) corrections are included 
that are necessary to achieve the theoretical precision for the total WW cross section of 
0.5% required for LEP2. In order to achieve a gauge-invariant description of the signal 
WW process and to employ the existing 0{a) electroweak corrections for the on-shell 
ly-pair production, we use the so-called leading-pole approximation (LPA). A double- 



pole variant of the LPA (usually referred to as DPA) was advocated in Ref. |]T9l and 
applied in a variety of calculations, see for example Refs. P,p!^-p!^ pO|] and Ref. for a 
review and further references. The important thing is that YFSWW3 and KoralW programs 
have a well established common part, which is the doubly ly-resonant (iyiy) process 
with the same universal radiative corrections. This allows us to combine the results of 
the two programs so as to achieve the desired theoretical precision for WW observables. 
Detailed descriptions of KoralW were published in Refs. [pT|-p^. Here, we present a 
documentation of the program YFSWW3; physical aspects together with some numerical 
results are discussed in Refs. p|-|6|,p0|]. Possible scenarios of combining the results of the 
two programs are presented in Refs. |l6|, pm|25|| . 

The outline of the paper is as follows. In Section 2, we define our notation and describe 
the physics contents of the program. In Section 3, we discuss the MC algorithm. In Section 
4, we present the structure of the program, important routines, etc. Details about the 
practical use of the program are given in Section 5. Finally, Section 6 summarizes the 
paper. Appendices contain useful technical information on the structure of the program, 
its input/output, etc. 
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2 Physics Contents of YFSWW3 



In this section we describe we define differential distributions used to calculate the cross 
section and generate events for the VT-pair production and decay process. These distri- 
butions we define quite completely, because they were not defined in every detail in the 
past publications. On the other hand we shall skip the detailed discussion of the physics 
models that they represent. This will be done in a separate publication First, we 



discuss our notation. Then we present the master formula and its ingredients. And finally, 
we give some details on the MC algorithms implemented in the program. 



2.1 Notation 




In YFSWW3 we adopt a notational convention similar to the one used in KoralW. The 
variables pi and p2 are the incoming e~ and four-momenta, respectively: 

Pl = (Ee, 0, 0, Pe), P2 = {Ee, 0, 0, -pe), (2) 

where is the beam energy and Pe = ^JE^ — m\ , with the electron mass mg. The 
variables Q\ and Q2 are the four-momenta of the W~ and bosons; their invariant 
masses are denoted by Mi and M2, respectively; gi, qi and mi,m2 are the four-momenta 
and masses of the W~ decay products; gs, and ms, m4 are the analogous four-momenta 
and masses of the decay products, with the first component of each pair corresponding 
to a fermion and the second to an antifermion. By fcj we denote the four-momenta of 
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radiative photons. See also fig. |I| for a pictorial representation. The components of a four- 
vector are denoted by p = {po,px,Py,Pz), with the metric convention of {1, —1, —1, — 1}. 
We also use the following Lorentz invariants: s = (pi +^2)^, Si = Ql = Mf and S2 = 

Ql = Ml 



2.2 Master Formula 

Here, we consider the process: 

e- + e^-^{W-^h + h] + {W^^h + h]+n^, (n = 0, 1, 2, . . . ), (3) 

where the radiative photons are emitted only at the VTM^-production stage. The photon 
radiation in W decays (FSR) in this version of YFSWW3 is done, in the 0{a'^) leading- 
logarithmic (LL) approximation, by the dedicated program PHOTOS - in the same way 
as in KoralW. Since this is done externally on fully generated events (it does not change 
event weights), we do not include the FSR in our discussion of the master formula. It 
should be remembered that the use of PHOTOS in YFSWW3 is a temporary solution before 
the advent of the YFS-exponentiated MC for W decays. 

As was already mentioned in the Introduction, we use the LPA to describe this process 
respecting the SU{2)l x U{1) gauge invariance. In YFSWW3, we implemented two options 
of the LPA, called LPAa and LPA^. The first one, which is the default (recommended) 
option, is based on the approach of Stuart where only the Lorentz scalar functions 



in the S-matrix are expanded about poles corresponding to unstable VT's. In the lowest 
order, the LPA^ matrix element coincides with the so-called CC03 (see e.g. Ref. |]15[ for 



an explanation of this notation) matrix element in the 't Hooft-Feynman gauge. In the 
following we shall use the notation "CC03" in that sense. The second option, intended for 



some dedicated tests, follows the method suggested in Ref. [15], where the whole S-matrix 



residuals are expanded about pole^; see Refs. for more details. In the following we 

shall describe the LPAq scheme explicitly; LPA;, can be obtained from this by taking the 
W^s on shell (i.e. Q\iQ\ — > M'^) everywhere, except for the Breit-Wigner denominators 
and the phase-space integration. 

In the framework of the YFS exclusive exponentiation (EEX), the cross section for 
the process (^ can be written as 

°° 1 /■ d?Qid^Q2 ^ d^qi " d^k- ( " \ 

n=0 ' Vl V2 Hi I \ i=\ ^ 

5^^^ (Ql - gi - g2) 5^^^ (Q2 - gs - q^) pT"\puP2; gi, • • • , g4, h,... , fc„), 

(4) 



■^We do not apply, however, the pole approximation to the phase space, as was done in the semi- 
analytical calculations of Ref. Il4| . 
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where p^'^^^ is the following fully exclusive multiphoton differential distribution: 

n 
i=l 

l + 5ciQi,Q2,Mw,Tw)] Yl PpB{{P,Q,q,k}'^'') 

pe{i,w}" 



y- ^Y'{pi,p2,Qi,Q2;k,) 



(5) 



where B is the contribution from the infrared-finite YFS residuals^ 
B{lp,Q,„,k}^')^[l+S^^f{[p,Q,„}^')] 

~t S {{p,Q,ki}'^p) 
t^j Sj{{p}^^, h) 2! ^ ~Si{{p}^^, k^) 



X ■ 



1_ f3i%K{{P,Q,(l}^'',ki,kj,ki) 



3! Sj{{p}^^, h)S,i{p}^^, fc,)5,(M^^ k 



where we use a short-hand notation of the type {p, Q, q} = {pi,P2, Qi, Q2, Qi, (I2, Qs, Qa}- 

The sum J2p Eq. runs over the photon partitions, i.e. all possible photon 
associations to the initial-state radiation (ISR) or the VT-state radiation (WSR). For n 
photons we identify each of the 2" partitions p with a vector p = {pi, . . . , where 
pi = I,W (/ and W denote the ISR and the WSR photons, respectively). The partition 
weight Pp is given by 



n Si{k,) n Sw{k,) n siik) n swik^) 

(7) 



=/ Pj=W pt=I pj=W 



EUSiiki) U Swikj) f[(Sj{k) + Sw{k 

where Si{k) and Sw{k) are the real photon infrared (IR) factors for the ISR and the 
WSR, respectively: 

a f Pi P2^ ^ 



Si{pi,P2,k) = - 

Att \kpi kp2 



OL ( Qi Q2 



SwiQi,Q2,k) = -^^^^^^ fcQ, 

It is worthwhile to stress that the partition dependence for the first two terms proportional 
to (3^^ and /Sj^'' cancels out, when truncated to 0{a) (as it should), thanks to the special 

•^In sums like we include all pairs of indices (i, j), so that i j , pi = I and pj = I. 
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normalization choice for the partition weights: ^ = 1. See also the definition 

pe{/,vy}" 

of P^'^ and the discussion on the partition-dependent "reduction procedure" TZ^ in the 
following. The purpose of the sum over partitions is to help the introduction of the LL 
corrections beyond the 0{a). 

In Eq. the two components proportional to jSj^^ and (3^^ explicitly represent the 
complete 0{a) contribution to the VrVT-production process, while the remaining /3's and 
A/?'s collect the second- and third-order ISR corrections in the LL approximation; see 
Section p.2.4| for more details. The dependence of the /3 functions on the photon partitions 



enters through the so-called reduction (extrapolation) procedures, which take into account 



the CMS energy shift induced by the ISR; for more details see Subsection |2.2.4 



2.2.1 YFS Infrared Factors 

The real photon IR factor 

includes the ISR, the WSR, and the interferences between them. 

The step function 6{k^ — k^), where = 61/5/2, cuts out the singular IR region, 
already included to all orders in the YFS form factor , where 

Y'{pup2, Qi, Q2; k,) = 2am'{pi,p2, Qu Q2) + 2aB{pi,p2, Qi, Q2] K). (10) 

The IR functions 3?-B' and B correspond to virtual and real photons, respectively, and 
are defined explicitly in Ref. |^. 

2.2.2 Coulomb Correction 

The correction 5c is the Coulomb effect, included both (two options) in its standard form 



as given in Ref. |28], and in terms of the so-called screened Coulomb ansatz of Ref. [29 
which is an efficient approximation of the non-factorizable corrections of Refs. [^D|-Q for 
singly-inclusive distributions. It reads 



5f'^'''\Qi.Q2.M^,T^) = ^ [1 - {1, (1 - /?)^} ^arctan (^^) 



:ii^ 



with 



/3 = ^1 - 2(si + S2)/SQ + [(Si - S2)/sq\\ P = ^(3^, SQ ={Ql + Q 




(12) 



4M, 



w 
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where Mw and Tw are the on-shell W mass and W width, respectively. The first term in 
the curly brackets in Eq. (pH]) corresponds to the standard (Std) correction, the second one 
to the screened (Scr) Coulomb correction. This correction is properly matched with the 
YFS form factor, where the Coulomb-like singularity appears in the virtual IR function 
After matching it with the full Coulomb correction 5ci the original YFS function 
3?-B is replaced by 3?-B', as described in Ref. 0. 



2.2.3 Anomalous Triple-Gauge Couplings 

The anomalous corrections to the triple-gauge couplings (TGC) WWV {V = Z,'j) are 
included multiplicatively through 

I A^AnTGC|2 

1 I XTGC _ K^'Born I /.oN 
^^"An — I . .SM 12 ' 
K^' Born I 

where A^Bom'^'^ Born-level matrix element with the anomalous TGCs, while A^lorn 
the corresponding SM matrix element. In YFSWW3 we use the same conventions (routines) 
for the TGCs as they are given in KoralW |2^. The anomalous couplings are included in 
the event generation only if an appropriate switch in the input parameters is on and the 
parameters corresponding to these couplings (in a given parameterization) are set up to 
non-SM values; see Appendix A for details. 



2.2.4 /3's for ly-Pair Production 

In the IR-finite /3 functions of Eq. (H), we see the superscripts 71^ in the four-momenta 
arguments like {p, Q, q}^^ or {p, Q, q, /c}^*". This means that these arguments are subject 
to the "reduction (extrapolation) procedures" . These reduction procedures are partition- 
dependent - this is why they are denoted with the subscripts p. The partition dependence 
enters through accounting for the ISR energy shift s s', where s' = (P1+P2 — X]p— / 
in the evaluation of the (3 functions. For more discussion on the meaning of the reduc- 
tion/extrapolation procedures, see Subsection Ol In YFSWW3, we use the reduction 



procedures of the YFS2 and YFS3 programs of KORALZ [p^-|36[l, with some modifications: 
in constructing the reduced W four-momenta, we take into account that, in general. 
Ml 7^ M2. We perform appropriate transformations of the decay products, that is 
we boost them to W^s rest frames before the reduction, and then make the inverse boost 
transformation to the reduced frame. In general, the effect (recoil) of the CMS energy 
shift s ^ s' is taken due to the ISR photons only. This procedure is well justified for 
the ISR corrections beyond 0{a) where we take {p, Q, q}^^ and for each photon we may 
say whether it belongs to the ISR or to the WSR. For the 0{a) specific corrections we 
cannot make a distinction, in the corresponding /5's, between the ISR and the WSR; 
the {p, Q, q, ki}^>^ in their arguments means that (i) the reduced {p, Q, ki}^^ are within 
the 3-body phase space and (ii) the photon ki is excluded from the energy shift s —>■ s' 



calculation. For more discussion on this non-trivial point see Subsections and 
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Let us now describe in detail the (9(a)-complete2 functions (3^^^ and (3[^\ The first 
one, which includes the 0{a) electroweak (EW) virtual corrections, reads as follows: 



Pi>\{p,Q}) = P'o'i{p^Q}) 



(0), 



l + 6i%{{p,Q};h)~Y{{p,Q};h 



(14) 



where is the lowest-order matrix element squared for the IV-pair production and 
decay; it is identical to the corresponding (3^^ function in KoralW |21,23|. In fact, the 
/3q^^ function comes in two variants. The first one, referred to as scheme (A) reads as 
follows: 



+ P^'\{p,Q}) 



1 + ^ISr(P1'P2) 

Q}; h) - Y{{p, Q}; K) - sSLipf^pf) 



(15) 



where ^jgl^ = (a / n)[\n(s / ml) — 1] is the LL pure QED 0{a) ISR correction, and the 
subscript means that the whole correction is calculated in the so-called scheme [0 
as explained in Ref. 0. Another option, referred to as scheme (B), similar to the one 
used in the program RACOON WW [|l^, is defined as follows: 



Piil{{p,Q})=P^'\{p,Q}) i + C(pi,P2) 



+ P^'\{p,Q}) 



5i'U{p,Q};K)-Yi{p,Q};K 



4^nipf,pf) 



a 



(16) 



where a is the fine structure constant and = {V^G sin^ 9w) / ■ Here, although 
the 0{a) corrections are calculated in the G^ scheme, their relative coefficient is a instead 
of ac ■ Our recommended option is scheme (A), which is inspired by the renormalization 
group equations (RGBs) P,P7|. 

The 0{a) virtual and soft photon corrections for the WW production stage are in- 
cluded in ^Ij+s- In the current version of the program, we use the calculations of Ref. |0 
for the on-shell pyVT-production. Since in our case the W^s are off-shell, we have to 
perform some transformations/extrapolations in order to employ these calculations. We 
do it in the following way. After performing the reduction procedure for the f3o function, 
as described above, we find cos6'efr and the velocity PeS of the W^s in the reduced frame, 
where 



/5eff = Vl-2(S1 + S2)/S'+[(S1-S2)/S']2, 
COS^eff = (pf'=-Qf^)/(b"fllQfl). 



(17) 



Then, for this value of (3eS we calculate the effective CMS energy squared for on-shell W's: 

4My(l-/5,2j. (18) 



ScS 



the early version of the present program caUed YFSWW2, corrections beyond the ISR were included 
only through the function Y' and 5- factors B . 
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Having {ses, PeS, cos9es), we construct the four-momenta (pf^ , pf^ , Qf^ , Qf^) for the on- 
shell process e+e~ W^W~ and use them to calculate the respective EW corrections. 
The same four-momenta are used to calculate the IR form factor Y according to Eq. (0); 
now the Coulomb-like singularity in the virtual IR function is kept, in order to cancel 
the corresponding singularity in S^Jg, and to calculate the 0{a) ISR correction: 

Opf,pf) = ^ (in ^-l). (19) 

Such a procedure allows us to use the EW corrections also in the VTiy-threshold region, 
however, in this region the LPA becomes less precise. Since there exists no other cal- 
culation with the 0{a) EW corrections for the VT-pair production and decay process in 
the lyVT-threshold region, we estimate, for the moment, the precision of YFSWW3 in that 
region at 2% conservatively, just as for the ISR approximation. This can be improved in 
the future, when more calculations/tests are done. 

Two comments are in order about the approximations used in the evaluation of 
Eqs. (|T§^): 



1. In the LPAa b, the one-loop corrections should be evaluated "on-pole", i.e. for the 
complex W mass: = —iMw^w- This would, however, require the analytic 
continuation of the on-shell one- loop corrections to the second Riemann sheet. We 
avoid this by using the approximation M ~ Mw in the LPA residuals. Since this is 
done for the 0{a) corrections, the error introduced by this approximation is of the 
order of {a/'7i)(T\Y/Mw) < 10~^, which is negligible for the precision we aim at. 

2. A coherent sum over W spins between the production and decay stages is included in 
Eq. ([T^) in the /5q°^ function, but not in the non- leading (NL) 0{a) EW correction. 
This should be sufficient for the LEP2 precision as the total EW correction in the 
LEP2 energy range is typically ~ 1-2%. However, the VT-spin effects can be included 
also in the NL 0{a) EW correction, if necessary. 

The 0{a) real photon ISR and photon emission in the lyiy-production stage are 
included through the following IR-finite function: 

/Si'^ {{p, Q, q, A:}^^ = D^^') {{p, Q, q, k}-^-) - S {{p, Q, q, fc}^'^) {{p, Q, g}^*') , (20) 

where D^^ is the differential cross section for the single hard photon radiation in the 
lyi^-production stage and the non-radiative M^-decays. For this we use the calculations 
of Ref. 0, where the W^s are taken on shell. In order to accommodate them in our off- 
shell process, we use general formulae for W^-polarization vectors in terms of their off-shell 
four-momenta in the spin amplitude calculations, and perform a coherent summation 
over the W's polarizations of the production and decay amplitudes. In fact, the hard- 
photon amplitudes of Ref. |^ are given in the massless-fermion approximation, so they 
become inaccurate when a photon is emitted close to a fermion direction. To correct 
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this, we supplement these amphtudes with the mass terms according to the CALKUL 
prescription of Ref. |^ . 



In Eqs. (ID and we use the IR-finite /? functions for the ISR /?q /gj^, /?} fgj^, P2 isr 

/^s^iSR' calculated up to 0(q;^) LL. They are exactly the same as in KoralW, see Ref. [p3| . 
We define in addition: 

A^(3) _^(3) _n{l) .ni3) _g{3) _g{l) .21) 

^A^OJSR ~ A^O.ISR A^OJSR' ^/^1,ISR ~ A^1,ISR A^1,ISR- y"^^) 

In Eqs. (§) and (|^), they are evaluated only for the ISR photons. In principle, these 
ISR residuals could be calculated for all photons - it is still consistent with the LL ISR 
approximation . 



2.2.5 Pretabulated EW Corrections 

We have implemented also the pretabulated (or "fast" ) version of the EW corrections. In 
this mode, a pretabulation of the NL 0{a) EW corrections 

^Ew(fe Q}) = SvUiP, Q}; K) - Y{{p, Q}- K) - 51sr(Pi,P2) (22) 

is done in the initialization stage of the program. Then, during the event generation, a 
linear interpolation of the pretabulated values is performed for each event. The pretab- 
ulation is done as a function of cos 6' of the W~ for fixed values of the CMS energy -y/s 
(set-up to the nominal energy) and of the W invariant masses (set-up to My/). This 
means, that in this approach, the ISR energy shift s — > s' as well as the off-shellness 
of the VT's is not taken into account (the latter prevents this method from being used 
below the VTH^-threshold) . In that respect, the pretabulated version of the EW correc- 
tions should be regarded as some approximation of the YFSWW3 "best" prediction. The 
main advantage of this version is that it is much faster in computing time than the exact 
( "best" ) one, because here the time-consuming EW library is invoked only at the begin- 
ning of the program run for some number of equally spaced points instead of by calling 
it for each generated event. We have checked that at -^i = 200 GeV, the results for 
the total cross section and main distributions of the pretabulated version differ from the 
corresponding results of the "best" one by only 0.1%, which is well below the required 
theoretical precision for the LEP2 data analysis. 



2.2.6 Fixed-Order Calculation 



By truncating the master formula (|[) at some power of the coupling constant a, one can 
obtain a fixed-order calculation, e.g. 0(q;°), 0{q^). 

At 0{q^) (the Born approximation), one gets the differential distribution 

pf{{p,Q.q})='PT{{P^Q^<l})- (23) 

At 0{q^), there are two contributions: (1) from the one-loop virtual corrections and 
the real soft photon radiation 

P^o^{{v.Q,q};K) = 'P^o\{v,Q.q}) \i + Cs{{v,Qy.K) \ , (24) 
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and (2) from the real hard photon radiation 

pS^)({p,Q,g},A;) = D«({p,Q,g},fc). (25) 

We have implemented the above fixed-order calculations in YFSWW3. They are available 
through the parallel weights, as explained in Table ^ of Appendix A. For 6^^^^ we use the 



■I. aiiu iui ij''^^ 

for more details. 



calculation of Ref. M, and for D\ the calculation of Ref. M; see also Subsection |2.2.4 



2.2.7 Radiative Corrections for W Decays 

Radiative corrections in W decays are not included in the above formulae explicitly. 
However, we include them in the overall normalization of the cross section by normalizing 
the W decay amplitudes to the W branching ratios (BR), which can be provided through 
the input parameters of the program. This means that in the (3 functions of Eq. (|^) for 
the {i,j)-th W~W^ decay channel, we use the normalization constant 

= -j^, (26) 

ref 

with 

_ awMw , . 

where 6w is the weak mixing angle, and aw = o:g^ in our recommended scheme, while 
we set aw = a in the so-called a scheme. Bi is the BR for the i-th W-decay channel, 
while Bref is the reference BR used for the normalization of the actual matrix element; see 



Ref. I^TI for more details. For the fully inclusive process (all channels), the normalization 
constant is 

Kii = (28) 

ref 

It turns out that the BRs calculated at the Born level in the scheme - the so- 
called improved Born approximation (IBA) - agree within 0.1% with the full one-loop 
results [^]. Such an IBA option is also included in the program - for this one needs only 
to supply (through the input parameters) the values of the CKM matrix elements. In 
that case 

{Bref '■ for leptons, 

3 |rCKM|2 (^^ ^ ^^iMM^ Bref ■■ for quarks, ^^^^ 

where V^^km ^^le CKM matrix element for the i-th 14^-decay channel, and as{M^) is 
the QCD coupling constant at the Mw scale. 

By construction, the real photon radiation in W decays, which is generated by PHOTOS, 
does not affect the overall normalization of the cross section - this is of course motivated 



by the KLN theorem [40,41 
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2.2.8 ISR Reference Differential Distribution 



Equations and (|^) define the so-called "best" physics model of YFSWW3; if the corre- 
sponding weight is used as a rejection weight in order to generate unweighted (constant- 
weight) events, then we say that YFSWW3 is run in the "best" mode. We want to combine 
the fully exclusive (multiphoton) differential distributions from YFSWW3 and KoralW, so 
that the resulting distributions contain both the LPA 0{a) EW corrections (missing 
in KoralW) and 4/-background diagram corrections (missing in YFSWW3), using event-per- 
event correction weights. In order to do this, we must define a certain simplified, auxiliary, 
fully exclusive, differential distribution da^, which is common to both programs, and im- 
plement the corresponding MC weight Wr in both programs. In YFSWW3 the weight Wr 
can be used also optionally, as a main rejection weight - this is the so-called "ISR mode" . 
Otherwise, wr is one of many alternative weights calculated simultaneously with the "best 
model weight"; it is available to the user, and/or it is used to calculate the correction 
weight that serves the purpose of introducing the missing 0{a) EW corrections in events 
generated by KoralW. 

In the reference distribution daji, only the numerically leading corrections, i.e. the 
ISR, the Coulomb correction, etc., are taken into account. This reference distribution 
includes the CC03 matrix element defined identically in YFSWW3 and KoralW. We also say, 
sometimes, that this is the distribution of the ISR-type model. The corresponding fully 
exclusive, reference, differential distribution, replacing p„ of Eq. (R), we define as follows: 



i=l 



X 



i=l 



Si {{p^, h 



i>j 



Sj{{p}^,h)Sj{{p}T^,k,) 



■y^ PsjSR {{P^Q^l}^^ki,kj,ki) 

,tr>l ^^({P}""' k^)SI{{p}'^, k,)Sj{{p}^, kj 



where we have made the following replacements with respect to Eq. (|^): 



Y' 

S 



Yi{pi,p2;e) 
SiiPi,P2,k), 



2^ 

71 



In^ - 1 lne + - In — - 1 + — 



(31) 



and we include the pure ISR LL /?'s only. The WSR is just absent from the above 
distribution, which means that there is only one photon partition (where all photons are 
associated to the ISR) with the weight Pp = 1. In the actual MC program, should be 
implemented with a MC weight in the operational mode when the WSR is switched off - 
through the parallel weights. Both of these possibilities are included in YFSWW3. However, 
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it can also be implemented as a MC weight when the WSR is switched on. In such a case, 
the sums over the /9's should extend to all photons. 

The above ISR reference cross section is also used in this and other works to define 
the so-called non-leading (NL) corrections 

c^NL = o- - o-R = y {da - dan) . (32) 

The subscript NL expresses the fact that these corrections are numerically smaller than 
the leading ISR corrections. It does not mean, however, that these corrections are free of 
logarithmic contributions, such as ln(s/M^), but these logarithms are numerically much 
smaller than the ISR basic logarithm ln(s/mg). 

As already indicated, the reference ISR differential distribution of Eq. (|30| ) (for the 
same CC03 matrix element) is also implemented in KoralW, and it has been checked 
numerically to a very high accuracy [§, ^ that its KoralW implementation agrees with 
that of the present YFSWW3 program. Having done that, it is possible to combine the 
results of the KoralW and YFSWW3 programs by reweighting the MC events, in order to 
include both the effects of the 4/ background and the LPA 0{q) NL corrections in the 
MC predictions for the H^-pair production and decay process. This is done in practice by 
reweighting the MC events produced by KoralW with the correction weight provided by 
YFSWW3. The correct ion to be used by KoralW is defined as 

^ ^i? ^ Pn'"*(PbP2;gi,--- ,q4,kl,... ,kn) 

where the numerator and denominators are defined in Eqs. (^) and (|30D. A few remarks 
are in due order. The above ratio is, in fact, calculated as the ratio of the corresponding 
MC weights. The YFSWW3 program has to be run in the mode with the WSR switched on, 
because p^"^^^ includes the summation over the ISR-WSR photon partitions, see Eq. (|^), 
through the multichannel reduction procedure TZp. Under these conditions the above 
correction includes the complete 0{a) correction, which can be used to correct the fully 
exclusive distributions generated by KoralW, on an event-per-event basis. As we shall see 
in Subsection [4.12|, the YFSWW3 program provides at present only a numerically efficient 



approximation of the above 5^] 



NL- 



2.3 YFS Reduction/Extrapolation Procedures 

It is characteristic and inherent feature of the YFS exponentiation that the IR-finite Pi 
distributions require extrapolation to a larger phase space with the additional "spec- 
tator" photons. For instance, the function /3o is in our case originally defined at the 
({^i; '?2, 93, 94}) = ({<?}) point of the 4-particle phase space and it has to be also defined 
at every ({g}, fci, . . . , kn) point of the (4 + ?T,)-particle phase space. We call the procedure 
of extending the domain of the (3o function /5o({g}) — *■ Po{{q}, ki, . . . , fc„) the "extrapola- 
tion procedure" . The only true limitation in the choice of the extrapolation procedure is 
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that /^(f ({g}, 0, . . . , 0) = Po{{q}) must hold. Typically, PoUq}) is defined in terms of dot- 
products of the four-momenta Pi-qj (or inner spinor products). Extrapolation may be done 
in a natural way by using exactly the same algebraic expressions in the larger phase space 
as were originally obtained from the Feynman diagrams in the smaller phase space. In the 
smaller phase space, the dot products obey certain relations due to the four-momentum 
conservation, which no longer hold true in the presence of the additional photons - this is 
the main mechanism in the realization of such an extrapolation. Another useful method 
of implementing the extrapolation (especially in the case when f3o{{q}) is implemented 
in the form of a black-box procedure, which requires as the input the four-momenta in 
the 4-particle phase space) is to project the {{q}, ki, . . . ,kn) point into the "reduced" 
point {q}^ using some kinematical manipulations (boosts, rotations, rescalings) on the 
four-momenta, and subsequently plugging it into /5(f ({g}, /ci, . . . ,kn) = f3o{{q}^). The 
above "reduction procedure" is, of course, less general than the extrapolation procedure. 
Last, but not least, let us note that the extrapolation procedure {{q}) {{q}, ki, . . . , fc„) 
is the reverse procedure to that of defining the residua at the IR-divergent points/poles 
({g}, ki, . . . , kn) — > ({g}) = ({g}, 0, . . . ,0) in the derivation of the resummation of the 
IR divergences to the infinite order (i.e. the derivation of the YFS exponentiation). The 
extrapolation procedure /3i({g}, kj) Pi{{q}, ki, . . . ,kj, . . . , kn) is quite analogous. It 
has to obey the relation /?i({g},%) /?f ({g}, 0, . . . ,0,/cj,0, ... ,0). Here all photons 
except the j-th one are the spectators. 

It should be stressed that the above freedom, due to a reduction/extrapolation proce- 
dure, is a well-known feature of the YFS exponentiation, already underlined in Ref. |^ 



and discussed in many works implementing the YFS exponentiation; see for example 



Refs. 0,0. The uncertainty of the results due to the reduction/extrapolation proce- 
dure in the C(a") YFS-exponentiated calculation is always at least of (9(a"+^); it does 
not infiuence or spoil the 0(0") perturbative "exact" contributions coming from the 
C(q;") Feynman- diagram calculations. On the contrary, a reasonable choice of the re- 
duction/extrapolation procedure may improve significantly the total precision, because 
the YFS exponentiation is then able to sum up efficiently higher-order effects beyond 
the "exact" 0{a"'), as was often seen |]4^-^. The above remark is also valid for the 
special "multichannel" variant of the reduction/extrapolation procedure described in the 
following. 

2.4 Multichannel YFS Reduction/Extrapolation 

It is useful sometimes to split the /3's into several components and introduce a different 
extrapolation/reduction procedure for each component. This is perfectly within the limits 
of our freedom. For instance in YFS3, in the centre of mass of the Z-boson, neither the 
initial beams nor the final fermion pairs have the momenta back-to-back. There are four 
possible definitions of the scattering angles in this frame. One defines f3o as a sum of four 
differential Born-level distributions, each for a different scattering angle. They might be 
summed up with the same weight 1/4, but, in fact, also a more sophisticated mixture is im- 
plemented 1^^. A similar multichannel extrapolation procedure is used in BHLUMI . 
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In the present work, we employ another kind of muhichannel extrapolation, introduced 



for the first time in the unpublished program BHWIDE ^6 . 

Let us elaborate on that method. Our exponentiation model is based on the real soft 
photon factor 5* and the corresponding virtual photon function to which photons 
emitted from the beams and from the W^s contribute coherently. We cannot, therefore, 
say whether a given photon is emitted from the initial state or from the W^s. We may, 
however, say something about it in a probabilistic way. If we split 5" as follows: 

S{k) = - A f ^ - ^ - ^ + = Si{k) + Swik) + SUk), 

An^ \kpi kp2 kQi kQ2 J ^^^^ 

Pi P2Y o /, ^ f Ql Q2^^ 



^^^^"^ \kpi A;p2/ ' ^^^^"^ \kQi kQ2 
we may define for each photon a probability 

= -^-^ (35) 

Si{k) + Sw{k) 

that it was emitted from the beams, and the probability pw{k) = 1 — pi{k) that it 
was emitted from the W^s. The above probabilities are a kind of LL concept, which 
has to be used with care. It may help to sum up higher-order QED corrections in the 
LL approximation. In particular, it may allow us to construct a multichannel type of 
reduction/extrapolation in which we may say for the spectator photons whether they were 
emitted by the beams (ISR) or by the W bosons (WSR). This is done as follows: for each 
phase-space point we split the differential distributions into 2" components, corresponding 
to 2" possible associations of n photons to the ISR or the WSR, the weight being the 
product P{k) of the corresponding pi{ki) or pw{ki). This product is exactly the partition 
weight in Eqs. (H) and (0). For each component we may apply a separate reduction 
procedure which "knows" whether a given spectator photon belongs to the ISR or the 
WSR. The above procedure is, in practice, simpler than is said above, because in the MC 
realization we start from the YFS3 differential distribution, in which Sint{k) is neglected, 
and the sum over the associations of photons to the ISR or the WSR is "randomized" . In 
other words, we generate only one of the 2" possible associations and therefore, for a given 
event, we work with only one component (association), for which we may immediately 
employ its "native" reduction procedure - that is we may use the information on the 
photon associations directly from the YFS3 MC generator; see also the next section for an 
additional discussion. Let us stress again that this powerful solution should be used with 
care, especially for the non-spectator photons in the (3 functions. 



3 Monte Carlo Algorithm 

In constructing the MC algorithm for the cross-section calculation and the event genera- 
tion, we start from the master formula of Eq. (^) and perform step-by-step simplifications 
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in the differential distribution, compensating them with appropriate weights. We do this 
until we reach a simple enough differential distribution to be generated using the stan- 
dard MC techniques. The MC algorithm of YFSWW3 for the Born-like WW production 
and decay is identical to the one used in KoralW for the so-called CC03 option ||21[], while 
the algorithm for the multiphoton radiation in the e~^e~ — > W^W~ process is based on 
the one implemented in the program YFS3 of KORALZ for the process e"*"e~ ff. 
Since W^s are much heavier than light fermions and, in general, have different invariant 
masses, we had to make a few modifications in the latter algorithm. The most important 
of them are: 

• generalization of the photon radiation formulae to allow for large and different 
masses of radiating particles^, 

• inclusion of interferences between the radiation from the initial state and that from 
the WW state, both in the YFS IR functions S and Y and in the IR-finite residuals 
p. 

The algorithm of YFS3, although used for a decade as part of other programs, was not doc- 



umented until recently. Its first full description can be found in two recent works [0,^ 
We therefore do not describe it here in detail, but refer the interested reader to the above 
two papers. We explain below only the main steps in the YFSWW3 algorithm that lead to the 
algorithms of YFS3 and KoralW, and then comment on the YFS reduction/extrapolation 
procedures. 

3.1 Main Steps in the Construction of the MC Weight 

The MC algorithm of YFSWW3 is based on the importance-sampling method. In con- 
structing it, we start from the master formula (H) and make a series of simplifications 
compensating them with appropriate weights until we reach a distribution that is simple 
enough to be generated using the basic MC methods (see e.g. Refs. [^Q), as we indi- 
cated above. Here, we describe only the main simplifications in Eq. (^) that lead to the 
formula that can be generated with the help of the MC algorithms used in the programs 
YFS3 and KoralW. They are: 

1. P functions: 

For the aggregate of the (3 functions we do the following replacement: 

B {{p, Q, g, A;}^-) ^ b, {{p, Q, g}^-) , (36) 
which is compensated by the "model" weight 

_ B{{p,Q,q,kr^) 



^This extension is also included in the YFS3 of the JCJC MC program Q 
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where bo = (3^^ is the Born-level differential cross section for the ly-pair production 
and decay process as given in Eq. (4) of Ref. |21[]. To evaluate the model weig ht w^ , 



we use the partition-dependent reduction procedures described in Subsection 2.2.4 



2. YFS form factor: 

We simplify this with the replacement: 

y'iPi,P2,Qi,Q2;K) — > Yj{pi,p2;e) + Yw{Qi,Q2;K), (38) 

which is compensated by the weight 

wy = exp[Y' - Yj - Yw], (39) 



where Yj is given in Eq. ( pT]) and IV is the YFS IR function corresponding to the 
photon radiation from the W~W~^ electric dipole: 

YwiQi, Q2; K) = 2am{Qu Q2) + 2aB{Q,, Q2, K). (40) 

The explicit formulae for the above and B are given in Ref. we do not 
repeat them here because they are rather lengthy. This simplification corresponds 
to neglecting the ISR-WSR interference terms in the YFS form factor. 

3. 5^ factors: 

In these we also neglect the ISR-WSR interferences: 

S{pi,p2,Qi,Q2,k) — > Si{pi,p2,k) + Sw{Qi,Q2,k), (41) 
and compensate this by the weight 



TT S{pi,p2,Ql,Q2,k) ^^2^ 

7=7 Si{pi,p2, k) + SwiQi, Q2, k) ' 



where Si and Sw are given in Eq. (|). 

After these simplifications, we obtain from Eq. (^) the formula for the first-level, 
"crude" cross section: 



If dQidQ2 dqi dk- r~ ~ ^ 

n=0 ^1 Ql [i^^ l^i ^ > 

5^'^ [vi +P2-Qi-Q2-f2 j ^^'^ -qi- q2) S^'^ {Q2 - gs - q^) 

pe{i,w}" 

(43) 
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The above distribution can be generated using the MC algorithms of YFS3 and KoralW. 
To make this more transparent, we shall write this formula in a slightly different but 
equivalent form. 

The product of the S-f actors in the square brackets of Eq. P3| ) together with the sum 
over the photon partition can be written as 



pe{i,w}" 



pe{/,vy}" 

n -^Sw{k,)e{k'l-h)bo{{p,Q,q}''^^) 



(44) 



where we have substituted Eq. for the partition weight Pp. 

After some algebra, and exploiting the Bose-Einstein symmetry for photons, the for- 
mula of Eq. (O) can be written as 



CCrude 



ni=Qnw=0-^ \ i=l / 

1 rl^k- 1 ^31. 



nw 



I - Qi - Q2 - E % 1 n ^ ^^'^ - ^1 - '^^'^ (^?2 - gs - 

V i=i / z=i ^' 



(45) 



where s' = P'^ is the CMS energy squared in the ISR-reduced (effective) frame, and nj 
and riw are the multiplicities of the ISR and WSR photons, respectively. 

To generate the ISR and WSR according to Eq. (^), we can use the MC algorithm of 
YFS3 (see Refs. |5^,^) for multiphoton radiation in the fermion-pair production, where 
the initial-final-state interferences were neglected. The only modification we had to do 
in this algorithm was to extend it to the case of heavy particles of unequal masses (see 
Ref. 10). Generating particular values of rij and nw corresponds to choosing a given 
photon partition for = n/ + nw in the sum of Eq. (0). This means that the sum 
over partitions is "randomized" in the program, i.e. not evaluated directly but generated 
using the MC methods. For the given (generated) photon partition, we then calculate the 
weights w^, wy and Wg. Such a solution can be regarded as a multibranch MC algorithm, 
see Ref. Q], where each branch corresponds to a particular photon partition. 



22 



3.2 MC Weights and Absolute Normalization 

In the process of simplifying the fully differential distribution of Eq. 

c^Best = ^ / dLips^iip; q}, h, . . . , K) Pn'^'^ip; q}, h,..., K) (46) 

n=0 ■ 

described in the previous section we have finished with the "crude differential distribution" 
^Crude j^gide the integral 



oo oo ^ ^ „ 

C^Crude= Y] Y] — T T / dUpS^{{p; q} , h, . . . pS^^({p;g},fcl,...,fcn), 



(47) 



n/=0 nw=0 



where n = n/ + ny/- Using the above notation, the total correcting weight is defined as 
follows: 

(1) Pn^'^*({P; Q'}; ^l! • • • ! ^n) p / .qx 

" PSS,({p;«}.fc.....W="'^'"'""^- 

In the actual program it comes in several versions, because for the purpose of the technical 
tests, and for the discussion of the physical precision, we need to switch off /on certain 
contributions in the full distribution p^'^'^*. 

However, in the actual MC program we do not actually generate the distribution 
^Crude^ but a Certain, more primitive, "primary" distribution pP''™. For its definition 
and full description, and the description of how events are generated according to this 
distribution, starting from the uniform random numbers, we refer the reader to the most 



complete documentation of the YFS3 Monte Carlo algorithm in Ref. [48|; Refs. |21] and [34 



can also be helpful. The weight correcting for the transition from the primary to crude 
distribution is given by 

^(.) ^ pZ:t({v:q}M.....K) 

and is fully defined in Ref. The total MC weight is of course equal to the product 

w = w^^^ w^'^^ (50) 

or 

w = w^°'^''^ (51) 

with 



^Model ^ ^Crudc ^ ^(2) ^-1^ (52) 



where w^°'^'^^ is called in the program the "model weight" (this weight and its many 
variants are provided in the WtSet array; see Section 5), and Crude jg Qg^Hg^j i^j^g "crude 
weight" (in the program it is equal to the product WtCrudl*WtCrud2; see Section 5). 
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The integrated cross section is given by the product of the average weight and the 
integrated primary cross section 



The integrated primary cross section we also call a "normalization cross section" , because 
it provides the absolute normalization for the whole MC calculation. 

4 Structure of the Program 

In this section we provide the reader with a brief guide of the YFSWW3 program. We shall 
describe its main routines, libraries and interfaces. We want to note here that the program, 
in the distribution version, is prepared for a UNIX/Linux- type operating system that 
supports directories and the make utility. The source code is distributed over a number 
of subdirectories, in order to make the structure of the program more transparent and 
easier to handle. Also a system of the Makefile's is provided, for compilation, execution 
and other auxiliary functions (e.g. clean-up). (Of course, this organization, which is in 
fact rather simple, can be avoided and the whole source code can be put into a single 
FORTRAN file.) In the following we give a review of all YFSWW3 subdirectories. We start 
from the directories that are most interesting from the user's point of view, then going to 
the ones with more technical contents. In the distribution package, all these subdirectories 
are located in the main YFSWW3 directory: yf sww3-l . 16-export. This directory also 
comprises two important files: README and RELEASE . NOTES, which we recommend the 
user look through before using the program. They contain some basic information about 
YFSWW3: how to compile/link and run the program, a brief documentation of the code, 
etc. 

4.1 demo — Demonstration Program 

The subdirectory demo contains a demonstration program in the file demo.f. Gener- 
ally, the user is supposed to provide his/her own main program; nevertheless, the file 
quoted here provides a simple example of such a program. It has a double role: (1) 
as a useful template, and (2) as a first cross-check that the MC generator YFSWW3 runs 
correctly on a given installation. The essential part of this program is a loop in which 
a series of MC events is generated. It also reads the input from a disk file, but no his- 
togramming is performed and most of the output comes from the generator itself. At 
the end of the program, a MC-integrated cross section of YFSWW3 is compared with the 
Born-level and ISR results from the semi-analytical program KorWan. The program is 
compiled/linked and executed, with the help of Makefile, for the input data set given 
in the files demo . input. The program also reads the default settings provided in the file 
data_DEFAULTS located in the directory data_f iles - some of these default settings are 
overridden by data from demo . input. The demonstration program is run for unweighted 
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events with the external hbraries PHOTOS, TAUOLA and JETSET switched on. The output, 
written into the disk file demo . output, can then be compared with the one provided in 
the file demo . output . linux. It was obtained on a PC Intel Pentium III under the Linux 
RedHat 6.1 operating system. 



4.2 data_f iles — Default Input Data 

The subdirectory data_f iles contains the input data file data_DEFAULTS, which gives the 
defaults settings for YFSWW3. It is identical with the one used in the program KoralW, so the 
two programs can be set-up from the same input data file! This considerably facilitates 
a combination the results of the two programs; see Ref. p5[ for more details. Some 



entries in this data file are dummy for YFSWW3 and are kept only for compatibility with 
KoralW, while some of them are specific to YFSWW3 and are treated as dummy parameters 
in KoralW. They are accompanied by appropriate comments. This file can also be used as 
a template to create the user's own input data file, where only those entries that will have 
different settings from the default data file may be included. The user input file should be 
read after the file data_DEFAULTS, so that the user settings can override the default ones; 
see the example in the demo . f file described in the previous subsection. This directory 
comprises also some older data files kept for backward compatibility. Their names contain, 
after a dot, the last two digits of the year of their creation. (The latest of these files is 
identical with the file data_DEFAULTS.) 



4.3 yf sww — Master Unit 

The subdirectory yf sww contains the actual Monte Carlo event generator. Subprograms 
that can be used directly are the following: 

• YFSWW_ReaDataX - the subprogram used to read, from the disk file, the default input 
data of YFSWW3 and subsequently the data of the user into the array xpar at the 
very beginning of the use of YFSWW3. 

• YFSWW_lnitialize - the subprogram that does all initializations of internal vari- 
ables. First, it sets-up the main parameters of the program (through the routine 
f ilexp) according to the values in the input data files, then it calls several initial- 
izers of the main internal modules of the generator, such as KarLud(-l, . . .), and 
of external packages, such as TAUOLA, etc. It prints out directly or indirectly all the 
input parameters. 

• YFSWW_Make - the most important subprogram of YFSWW3. It generates single MC 
events. Functionally, it is a high-level management subprogram in the event gener- 
ation. It invokes (through the routine yf sww3) other routines that perform specific 
tasks, such as: generation of ISR photons (KarLud), generation of radiative photons 
from W^s (KarFin), generation of the W^s and final 4/ four-momenta (WW_Presam), 
evaluation of the Coulomb correction (WTCoul), calculation of the matrix element 
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(Model), of the YFS form factor, of the interference corrections, etc. It builds up 
the total MC weight from the weights supplied by these routines. This weight is 
then returned, in the weighted-event mode, as the main event weight, or used in the 
rejection loop (in the routine yf sww3), in the unweighted-event mode, for construct- 
ing the weight = 1 event. It also calls PHOTOS, which generates photon radiation 
in the W decays, TAUOLA for decaying of the final-state r's (if they appear), and 
JETSET to perform fragmentation/hadronization of the final-state quarks. It keeps 
track, through special monitoring routines, of all the information necessary for the 
final results from the MC event generation. 

• YFSWW_Finalize does all final bookkeeping, including the calculation of the inte- 
grated (total) cross section. It prints a summary output for the whole sample of 
generated MC events. 

The first subprogram is located in the file readata.f , while the following three are in the 
file yf sww3 . f . 

This directory contains also a collection of utility routines, which can be used, for 
instance, to access information on particle flavours and four-momenta, on event weights, 
etc. They are located in the file ww_get . f (see the comments in this file for more details). 

4.4 model — Matrix Elements 

The routines for the "model" weights calculations are located in the subdirectory model. 
The master routine here is the subroutine Model (in the file model. f); by calling some 
other routines, this evaluates the "model" weights w^°'^'^^ as defined in Eq. (|52|). 

1. For the Born-level matrix element we use the same routines as in the program 
KoralW - they are collected in the file born.f . The same again are used for the ISR 
corrections, located in the file betas. f. 

2. The 0{a) virtual and real soft photon corrections are evaluated with the 
help of a special interface - the routine VirSof - in the file virsof .f . In the first 
call, this routine makes all the necessary initializations, and then calculates the EW 
corrections by calling appropriate routines from the EW library for each generated 
event. The pretabulated (or "fast") version of the EW correction is invoked through 
the routine DnlFast (in the file dnlf ast .f ). This routine makes a pretabulation of 
the EW corrections in the initialization stage of the program and then, during the 
event generation, it performs the linear interpolation of the pretabulated values to 
calculate the EW corrections for each event as described in Section p.2.5| . 

3. The 0{a) hard photon matrix element in the WW production, which con- 
tributes to the function Pi, see Eq. (PP]), is calculated with the help of the routine 
eewwg, located in the file eewwg.f. This routines evaluates appropriate spin ampli- 
tudes for the process e~^e~ — > W^W~'y according to the formulas of Ref. 0, which 
are then combined coherently with the spin amplitudes for W decays. 
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In order to calculate all these contributions to the (3 functions, the subroutine Model calls 
the routines that perform appropriate "reduction" procedures as described in Section 2. 
Except for calculating the contributions to the YFS exponentiated cross section, this 
routine provides also the weights corresponding to the strict 0{a^) and C(a^) calculations. 
In the end, the subroutine Model fills the array of weights with the weights corresponding 
to: (i) several variants of the SM models (for example different orders of QED ISR) and 
(ii) various components in the differential cross section for a given variant of the SM (for 
example contributions from various /?'s). The subroutine Model returns the best "model" 
weight. 

The routine WTCoul - in the file coulco.f - for the Coulomb (both "standard" and 
"screened" ) correction calculation is also located in this directory. 

4.5 ewe — Electroweak Library 

The subdirectory ewe contains the routines for the calculation of the electroweak (virtual 
plus soft photon) corrections in the on shell WW production process according to Refs. |^, 
. For the evaluation of one- loop integrals they use the package FF of Ref . |^ , which 
is located in a separate directory (see below). 

4.6 interf s — Interfaces to External Libraries 

Interfaces to some external libraries are collected in the directory interf s. The most 
important routines are: 

• inietc - in the file tauola_photos_ini . f - sets up some parameters of PHOTOS and 
TAUOLA. 

• tohep - in the file hepf ace . f - fills in the /HEPEVT/ COMMON block, calls PHOTOS for 
photon radiation in W decays and TAUOLA for r decays. 

• tohad - in the file lundf ace . f - calls JETSET for the quark fragmentation/hadronization. 

Note that the COMMON block /HEPEVT/ is expected to contain single-precision (REAL*4) 
variables. 

4.7 wdeclib — External Libraries for W Decay 

The subdirectory wdeclib provides the following programs: 

• PHOTOS 1^ for photon radiation in the W decays (up to two photons), 



TAUOLA 1 52 1 for r decays with radiative corrections. 



JETSET version 7.4 [|53[ for quark fragmentation/hadronization, 
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located in the files: photos. f, tauola.f and jetset74.f, respectively. 

In YFSWW3 1 . 16 we use the recently modified PHOTOS, where also the photon radiation 
from quarks can be activated with the help of a special input parameter switch; see the 
tables with the input parameters in Appendix A. Note that in the distribution version of 
TAUOLA the parameters in the r-decay modes are not adjusted to the recent experimental 
data. We recommend the user to replace this version of TAUOLA with the one accepted 
within his/her own collaboration. The technical update described in Refs. p4|,|55| will 
resolve this inconvenience in the future releases of the program. 



4.8 glib — Histogramming and Numerical Library 

A handy FORTRAN histogramming package, GlibK ^6|], is provided in the subdirectory 
glib (the file glibk.f). It is used by YFSWW3 both for hard-coded internal bookkeep- 
ing and for some optional "external" tests. The package is similar in its usage to the 
classic HBOOK of CERNLIB. This subdirectory also contains the useful numerical package 
yf slib.f with the routines for weight monitoring, random-number generation, numerical 
integration, one- and two-dimensional adaptive MC sampling, Lorentz transformations, 
etc. 



4.9 semian — Semi- Analytical Program KorWan 

The subdirectory semian contains the package KorWan for the semi-analytical calculations 
of the cross section for the off-shell WW production and decay at the Born level and with 
the ISR corrections (using the structure- function formalism). Also the Coulomb effect 
is included in KorWan in both versions: the "standard" and the "screened" ones. The 
package is described in detail in Refs. ETl-pS. It can be used for a quick consistency 
check of the YFSWW3 results. 



4.10 ff - FF Package 



In the subdirectory f f is located the FF package |51| for the evaluation of one- loop integrals 
needed to calculate the 0{a) EW virtual corrections. 



4.11 dok — Related Papers 

In the subdirectory dok, we include some of our papers related to YFSWW3. They are in 
the form of gzip'ed PostScript files. 



4.12 rewt — Reweighting Tools 

The subdirectory rewt contains routines that are not needed for the standard MC event 
generation but for various event reweighting purposes. We provide here the tools for 
two kinds of reweighting: the reweighting (1) of the KoralW generated events to correct 
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for 0{a) NL effects in WW production, and (2) of the YFSWW3 generated events to take 
into account effects due to some input parameters (typicaly H^-mass) changes. All these 
tools are described in detail in the following; see also the file README provided in this 
subdirectory. 



4.12.1 Reweighting KoralW Events by YFSWW3 

The main tool for this kind of reweighting is the function YFSWW_WtNL, located in the 
file rewt_K.f. For the events, in terms of flavours and four- momenta of the final-state 
fermions and the ISR photons, provided through its parameters, this function calculates 
and returns the corresponding correction weight from YFSWW3. This weight is a ratio of 
the fully differential distribution. In practice, we calculate the ratio of the MC weights 
corresponding to two distributions - this is more convenient. We can do it because the 
common crude distribution entering both weights cancels out, for a given number of real 
photons of Eq. (^, to the reference distribution in the "ISR approximation" defined in 
Eq. (^), which is the same as in KoralW in the CC03 mode (Eq. (4) of Ref. [^). It is 



defined in Eq. (|33|). In the reweighting mode of YFSWW3, we implement an approximate 
version of this weight: 

_ pT'-{{p,Q,qr',h,...,K) _ 

where Pn^^^~ is a variant of p^^^^ in which the sum over the photon partitions is reduced to 
only one term, in which all photons are associated to the ISR with the partition probability 
= 1, and TZj means the respective reduction procedure. Unfortunately, we loose in 
this way a (small) part of the genuine 0{a) corrections. The full 0{a) correction can be 
included only if YFSWW3 is run in the event-generation mode, and the correction weight 
for reweighting (accounting for the 4/ background) is supplied from KoralW. (This leads, 
however, to large MC weight fluctuations, see below.) We have checked, however, that 
the bias introduced in the LEP2 energy range by this approximation is within 0.1%, both 
for the total cross section and for distributions. This can be explained by the fact that 
the W^s, as heavy and slowly moving particles, do not radiate much at these energies. Of 
course, the main weight of YFSWW3 based on p^^*^* does not have this bias. In this way, by 
assigning this correction weight to the KoralW events, one can include, in a simple way, 
the 0{a) NL correction into the KoralW results. How to do it in practice, we explain 



in detail in Ref. 25 



4.12.2 Reweighting YFSWW3 Events by YFSWW3 

This kind of reweighting can be useful in some data analyzes, e.g. when one has accu- 
mulated a large sample of simulated events for a given input parameter set-up and one 
wants to assess the effects of changing some of the input parameters, but without hav- 
ing to repeat the full simulation with the new input parameters. The tools for such a 
reweighting are collected in the file rewt_Y.f. We provide two methods: "exact" and 
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"approximate" . In the exact method the partition vector for photons is stored along with 
their four-momenta, while in the approximate method the partition information is not 
recorded. 

(a) The exact method: 

The tool for the exact method is SUBROUTINE YFSWW_WME(WtME) , where WtME corresponds 
to the weight w^°'^'^^ of Eq. (0). It requires storing on the disk/tape the entire content 
of four COMMON blocks: /MOMINI/ , /MOMFIN/ , /MOMDEC/ and /DeChan/ during the event- 
generation run of YFSWW3. We assume that the appropriate tools for writing the above 
information on the disk/tape as well as for reading it in the reweighting run of the program 
are provided by the user of the code (see the source code for more details). 

How to use it: 

1. The primary-run generating/storing events: 

• Run YFSWW3 and store on the disk/tape the whole content of the COMMON blocks: 
/MOMINI/, /MOMFIN/, /MOMDEC/, /DeChan/. 

This is done for some input parameters set-up (A). 

2. The secondary-run reweighting events: 

• Initialize the program with the input parameters set-up (A) used for the event 
generation. 

• Instead of generating events, call YFSWW_WME(WtME_A) for all stored events and 
save the computed matrix element weight WtME_A (in some arbitrary units). 

• Initialize the program with a different input parameters set-up (B). 

• Calculate the weight WtME_B and save as above. 

• For each event, calculate the ratio Wt = WtME_B/WtME_A and use it as the weight 
for reweighting the original sample of events. 

• This procedure can be repeated for any new input parameters set-up (B). 

Example of running the program in the reweighting mode: 

CALL YFSWW_ReaDataX(' ./data_DEFAULTS' ,1,10000, xpar) ! reading general defaults 
CALL YFSWW_ReaDataX(' ./user. input' ,0, 10000, xpar) ! reading user input 
CALL YFSWW_Initialize(xpar) ! initialize generator 

NevTot = 1000 ! number of events to be reweighted 

DO lev = 1, NevTot 

CALL YFSWW_WME(WtME) ! calculate matrix element weight for stored events 
ENDDO 

For more information about running YFSWW3, see Section ^. 

(b) The approximate (simpler) version: 

In the case when one does not have all the information needed for the exact (recommended) 
method of reweighting (because of the lack of disk/tape space or some other reasons), one 
can use 
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WtME = YFSWW_WME_Smpl(Iflav,pf l,pf2,pf3,pf4,Phot,Nphot) 

instead of CALL YFSWW_WME(WtME) . In this case one needs to store for each event only 
the flavours and four-momenta of the final-state fermions and the number and the four- 
momenta of all photons. Important: They should come from YFSWW3 itself, not from 
PHOTOS after the FSR. These four-momenta and flavours are available from two COMMON 
blocks: 

COMMON / DECAYS / If lav(4) ,amdec(4) 

COMMON / MOMDEC / pf 1 (4) ,pf 2 (4) ,pf 3(4) ,pf 4(4) ,Phot (100 ,4) ,Nphot 

(the fermion mass array amdec(4) does not need to be stored). They can be alternatively 
accessed with the help of two getter-routines: 

SUBROUTINE YFSWW_Get4f (f lav , pi , p2 , p3 , p4) 
SUBROUTINE YFSWW_GetPhot All (NphAll , PhoAll) 

located in the file ww_get . f in the subdirectory yf sww. The first one provides the flavours 
and four-momenta of the final-state 4-fermions, while the second one provides the number 
and four-momenta of the radiative photons (excluding the ones from PHOTOS). 

The FUNCTION YFSWW_WME_Smpl fills the appropriate COMMON blocks and calculates (re- 
turns) the matrix element weight (corresponding to WtME in the previous method). This 
method is approximate because it does not use the information on the photon partitions 
but assigns all the photons to the ISR in performing the partition-dependent reduction 
procedure (described in Section ^.2] ). Our tests show, however, that the differences be- 
tween this and the exact method are numerically negligible, e.g. < 0.01% for the total 
cross section at y/s = 200 GeV. 



4.12.3 Reweighting YFSWW3 Events by KoralW 

Similar reweighting tools are provided in the new version of KoralW They can be used 
for reweighting events generated by YFSWW3 with the weight provided by KoralW in order 
to correct for the missing 4/-background contribution to the signal WW process. The 
relevant reweighting routines in KoralW require as an input the events from YFSWW3, which 
are generated according to the ISR-type reference differential distribution. Such events are 
provided by YFSWW3 in the standard generation mode, with the help of the getter-routine 
YFSWW_GetEvtISR, located in the file ww_get.f of the directory yfsww. The correction 
weight provided by KoralW can be used to reweight the constant-weight (unweighted) 
events from YFSWW3 or to correct the YFSWW3 weights in the variable-weight (weighted- 
event) mode; see Ref. pS] for more details. 



5 How to use the Program 

In this section we provide a short guide of how to use the current version of YFSWW3, and 
we describe the input parameters of the program and its output. 
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5.1 Principal Entries of YFSWW3 

The principal entries of the YFSWW3 package, which the user has to call in his/her applica- 
tion in order to generate a series of MC events, were already listed and described briefly 
in Subsection W^. Here we shall add more information on their functionality. The calling 



sequence constituting a typical Monte Carlo run will look as follows: 

CALL YFSWW_ReaDataX(' ./data_DEFAULTS' ,1,10000, xpar) ! reading general defaults 

CALL YFSWW_ReaDataX(' ./user. input' ,0, 10000, xpar) ! reading user input 

CALL YFSWW_Initialize(xpar) ! initialize generator 

DO loop = 1, 10000 ! loop over MC events 

CALL YFSWW_Make ! generate single MC event 
EMDDO 

CALL YFSWW_Finalize ! final bookkeeping, print 

CALL YFSWW_GetXSecMC(XSecMC,XErrMC) ! get total cross section 

In the first call of YFSWW_ReaDataX, default data are read into the REAL*8 array xpar (10000) . 
The YFSWW3 has almost no data hidden in the source code (this is not true for TAUOLA 
and JETSET). The file data_DEFAULTS, which is read first, is placed in the distribution 
subdirectory data_f iles. This file provides necessary initial default values of all input 
parameters. The user should never modify it. It can be copied to a local directory or, bet- 
ter, a symbolic link should be created to the original data_f iles/data_DEFAULTS. This 
file is quite sizeable and the user is usually interested only in changing some subset of 
these data. In the second call of YFSWW_ReaDataX, the user can overwrite the default data 
with his/her own smaller set of input data, which are placed in the user, input file. For 
example, the simplest input data, defining only the CMS energy, would look like this: 

BeginX 

*<ia>< data >< comments > 

1 190d0 CMSEne =CMS total energy [GeV] 

EndX 

As we see, data cards begin with the keyword BeginX and end with the keyword EndX. The 
comment lines are allowed - they begin with * in the first column. The data themselves 
are in a fixed format, with the index i of the array xpar (i) followed by the data value and 
a trailing comment. The example of the input data set for the demonstration program 
demo . f in the subdirectory demo provides a useful template for typical user's data. The 
complete set of data in data_DEFAULTS is described in detail in Tables see Appendix 
A. Obviously, the user is interested in manipulating only some of them and will retain 
the default values in most of the cases. 

The YFSWW_Initialize is invoked to initialize the generator. It reads the input data 
from the array xpar, prints them out and sends them down to the various modules and 
auxiliary libraries. The programs have to be called strictly in the same order as in the 
above example. At this point one is ready to generate the series of MC events. The 
generation of a single event is done with the help of YFSWW_Make. After the generation 
loop is completed, we may invoke YFSWW_Finalize, which does the final bookkeeping. 
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prints out various pieces of information on the MC run, and calculates the total MC 
integrated cross section (in pb). In order to obtain this cross section the user may call 
the routine YFSWW_GetXSecMC(XSecMC,XErrMC). 



5.2 Input/Output Parameters 

As we explained in the previous section, the input parameters enter through the xpar 
array, being a parameter of the routine YFSWW_Initialize. Their meaning is explained 
in Tables of Appendix A. 

The principal output of YFSWW3 is the Monte Carlo event, which is just a list of 
final-state four-momenta (in GeV) and flavours, encoded in the standard /HEPEVT/ event 
record. In the present version, we still provide the REAL*4 version of /HEPEVT/ of dimen- 
sion 2000. If the user is interested in the parton momenta before hadronization, then, in 
addition to /HEPEVT/, they are available (see also Table ||) through 



CALL YFSWW_Get4f (pl,p2,p3,p4) 

CALL YFSWW_GetBeams(ql,q2) 

CALL YFSWW_GetPhotAll (NphAll , PhoAll) 



get final 4f flavours and 4-momenta 
get beam 4-momenta 

get photon multiplicity and 4-momenta 



Alternatively, all the four-momenta are available from the internal COMMON blocks /MOMWWP/ 
and /MOMDEC/, see Tables | and | for details. 

In the case of a MC run with weighted events, the user is provided with the main 
weight (see also Table ^) through 

CALL YFSWWW_GetWtMain(WtMain) ! get main Monte Carlo weight 

Of course, for unweighted events WtMain=l. For special purposes the user may also be 



interested in auxiliary weights, which are provided (see also Table |T0]) with the help of 



CALL YFSWW_GetWtAll(WtMain,WtCrud,WtSetAll) ! get all Monte Carlo weights 

where REAL*8 WtSetALL(lOO) is an array of the weights described in Table ^ The 
definition of the main weight WtMain depends on the input switches KeyCor and KeyLPA, 
see Table ^; for instance, the best model for KeyCor = 5 and KeyLPA = corresponding 
to Eq. (H) is defined as follows: 

WtMod = Wtcrul*Wtcru2*(WtSetAll(41)-WtSetAll(2)+WtSetAll(4)) . 

Alternatively, the auxiliary weights WtSet defined in Eq. (^) are available from the 
COMMON /WGTALL/, see Table |. The total auxiliary weight should be defined as: 

WtAuxEWC(i) = WtCrul*WtCru2*WtSetAll(i) 

for 0{a) EW corrections, and 

WtAuxISR(i) = WtCrul*WtCru2*WtSetAll(i) 

for the ISR corrections. The corresponding integrated cross section is simply obtained by 
multiplying the average of the total weight by the normalization cross section of Eq. 



(or the "primary cross section" in the terminology of Ref. which is provided through 
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CALL YFSWW_GetXSecNR(XSecNR,XErrNR) ! get normalization x-section 

Note that we expect the user to exploit the auxihary weights only in the variable-weight 
operation mode. It is, however, not impossible to use them also for the constant-weight 
events. In such a case we recommend the user to contact the authors for more instructions 
on how to do it correctly. 

The complete description of the post-generation output parameters from YFSWW_Finalize 
is collected in Table |TI|. 

5.3 Printouts of the Program 

In this section we describe a printout of the demonstration program demo in the "best" 
mode, shown in Appendix B. The printout starts with the detailed specification of the 
actually used input parameters. Also, logos of some of the activated subprograms (KarLud, 
KarFin) and libraries (PHOTOS, TAUOLA, JETSET) are printed here. Next, the printout of 
one full event (in the standard PDG convention) is shown (in the actual demo . output file, 
five events can be seen). Then, the printouts from the post-generation mode, i.e. CALL 
YFSWW_Finalize, appear. First, one can see summary reports from the subprograms 
KarLud and KarFin. Their meaning is rather technical, so the average user of YFSWW3 
does not have to worry about them (unless something unusual appears there). One thing 
that the user may check from time to time there is the entry B5 in the KarFin window 
B: the average value of the WCTRL weight should be equal to 1 within the statistical error 
(if it is not, please notify the authors!). After these technical printouts from the YFSWW3 
subprograms, there are several summary reports from TAUOLA on the r decays. The final 
reports of the main YFSWW3 MC generator are collected in three windows: A, B and C. 

The window A contains the most important information from the user's point of view. 
It provides the values of: 

AO : the CMS energy in GeV; 

Al : the best-order total cross section with its absolute statistical error for the generated 
statistics sample (in pb) ; 

A2 : the relative error of the above cross section; 

A3 : the total number of generated events; 

A4 : the number of accepted events; 

A5 : the number of events with negative weights; 

A6 : the relative contribution to the total cross section from the events with the negative 
weights; 

A7 : the number of overweighted events; 

AS : the relative contribution to the total cross section from the over-weighted events; 
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A9 : the value of the maximum weight for event rejection; 
A 10 : the average weight. 

The numbers in the entries A5,A6,A7,A8 should be as small as possible. In particular, the 
relative contributions to the cross section from the negative weights and the overweights 
should be much smaller than the expected accuracy of the MC calculations. If there is a 
large contribution from the overweighted events, one should try to increase the value of the 
maximum weight WtMax, see Table |l|. In the case of a large negative weight contribution, 
one should contact the authors. The ratio of the average weight (entry AlO) to the 
maximum weight (entry A9) shows the efficiency of the MC algorithm in the unweighed 
events generation, i.e. the event acceptance rate. In the example given in the demo 
program, it is ~ 20%. 

The window B is devoted to the technical information on the ISR corrections in different 
orders in a within the YFS framework. This information is quite important, because it 
shows how big the contributions of subsequent orders of the perturbative series are, and 
thus allows us to estimate the missing higher-order effects. The entries B1-B4 provide 
the values of the total cross section (in pb) in the orders 0(Q;°)exp^C'(a'^)exp of the ISR, 
while the entries B15-B17 give the differences of these values in subsequent orders. The 
entries B5-B14 contain the information on the contributions from the individual YFS ISR 
residuals PijsR^ also in different orders in a. The differences between various residuals 
and various orders for a given residual are provided in the entries B18-B26. 

The window C contains the information on the exact (9(a^)exp corrections to the WW 
production stage. The first entry (Cl) shows the total C(a^)exp cross section for the exact 
EW corrections calculation, the second one (C2) shows the same for the approximate 
("fast") EW corrections, and the third one (C3) contains the 0(a^)cxp LL cross section. 
The second part of this window shows the differences between the exact and approximate 
EW corrections (C4), and between the exact and the LL results (C5). The final part 
(entries C6-C8) contains the differences of various contributions. All these values are 
given in pb. 

This completes the description of the output of YFSWW3. The remaining entries shown 
in the demo . output file (e.g. histograms) are produced by the demo main program. 

6 Summary 

We presented the Monte Carlo event generator YFSWW3 version 1.16 for the combined Im- 
pair production and decay process. It includes the complete 0{a) EW corrections in the 
WW production process in the leading-pole approximation, in addition to other numeri- 
cally sizeable but physically less interesting effects, such as QED ISR effects up to 0{a^) 
in the LL approximation. The program is precise enough to obtain a Standard Model pre- 
diction for the total cross section and any distribution at LEP2, understanding that the 
contribution from the background diagrams is included in the calculation with the help of 
KoralW (or removed from the data). YFSWW3 includes programming tools to communicate 
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with KoralW in the process of generating a single MC event. The YFSWW3 program was 
also tested and is applicable for the energy range up to 1.5 TeV (the LC/TESLA range). 
The next possible improvement in the program is the introduction of the complete YFS 
exponentiation in the MC simulation of the W decays. 
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A Program Parameters and Their Settings 



Parameter 


Position and meaning 


CMSEne 


xpar (1) 


C=200.0): y/s, centre-of-mass (CMS) energy [GeV] 
C=l. 16639d-5): Gp, Fermi constant [GeV~^] 


Gmu 


xpar (2) 


alfWin 


xpar (3) 


[=128 . 07d.O) : I/ck^k inverse QED coupling constant at My^ scale 


aMaZ 


xpar (4) 


[=91 . 1882) : Mz, mass of Z boson, [GeV] 


GammZ 


xpar (5) 


[=2.4952): Tz, width of Z boson [GeV] 


aMaW 


xpar (6) 


[=80.419): Mw, mass of W boson [GeV] 


GammW 


xpar (7) 


[=-2.120): r^^, width of W boson [GeV], for (gammw < 0) Tw 

is recalculated from G^, M-^r and as 


VVmin 


xpar (8) 


[=ld-6): Minimum ?;-variable (dimensionless) , IR cut-off 


VVmax 


xpar (9) 


[=0 . 99) : Maximum value of v-variable 


WtMax 


xpar (10) 


(=2.0): Maximum weight for rejection, for wtmax < redefined 

inside the program 


aMH 


xpar (11) 


(=115.0): Higgs mass [GeV] 


aGH 


xpar (12) 


(=1.0): Higgs width [GeV] 


alpha_s 


xpar (13) 


(=0.1185): QCD coupling constant 




xpar (14- 19): Dummy parameters in YFSWW3 (specific to KoralW) 



Table 1: The list of input parameters of the YFSWW3 generator in the xpar vector. The default 
values are in brackets. 



PfirauKM cr 


Pt)sitiou aud mcvuiing 


KeyCor 


xpar (2001) (=5): Radiative corrections switch 




=0 


Born 




=1 


Above + ISR 




=2 


Above + Coulomb correction 




=3 


Above + Full YFS form factor for WW production 




=4 


Above + Radiation from WW 




=5 


Above + Exact 0{a) EW corrections in WW production (BEST!) 




=6 


Same as 5 but with pretabulated EWC (approximate but faster!) 


KeyLPA 


xpar (2002) (=0) : LPA mode switch 

=0/1: LPAa/LPA(, (LPA^ recommended, LPA^ for test only) 



Table 2: The list of input parameters in the xpar vector specific to the YFSWW3 generator 
(not used in KoralW). The default values are in brackets. 
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Parameter 


Position and meaning 




xpar(21-57): Values of TGC's 


: Set 1, most general set - complex numbers. 




default values are wild random, 


not shown 


gl(l) 


=DCMPLX(xpar(21) ,xpar(31)) 


= , for WWZ vertex 


kap(l) 


=DCMPLX(xpar(22) ,xpar(32)) 


= for WWZ vertex 


lam(l) 


=DCMPLX(xpar(23) ,xpar(33)) 


= Xz, for WWZ vertex 


g4(l) 


=DCMPLX (xpar (24) , xpar (34) ) 


= gl, for WWZ vertex 


g5(l) 


=DCMPLX(xpar(25) .xpar (35)) 


= g^, for WWZ vertex 


kapt(l) 


=DCMPLX(xpar(26) , xpar (36)) 


= kz, for WWZ vertex 


lamt(l) 


=DCMPLX(xpar(27) , xpar (37)) 


= Xz, for WWZ vertex 


gl(2) 


=DCMPLX(xpar(41) ,xpar(51)) 


= gf, for WW 'J vertex 


kap(2) 


=DCMPLX(xpar(42) , xpar (52)) 


= Kg, for WW J vertex 


lam(2) 


=DCMPLX(xpar(43) , xpar (53)) 


= Xg, for WW^ vertex 


g4(2) 


=DCMPLX(xpar(44) , xpar (54)) 


= g^, for WW^ vertex 


g5(2) 


=DCMPLX(xpar(45) ,xpar(55)) 


= 55 , for WW 'J vertex 


kapt(2) 


=DCMPLX(xpar(46) , xpar (56)) 


= kg, for WW^ vertex 


laint(2) 


=DCMPLX(xpar(47) , xpar (57)) 


= Xg, for WW^ vertex 




xpar (61-65): Values of TGCs: Set 2, see CERN 96-01, Vol. 1, p. 525 


delta_Z 


=xpar(61) = 8z 




X .gamma 


=xpar (62) = 




x_Z 


=xpar (63) = xz 




y gamma 


=xpar(64) = 




y-z 


=xpar(65) = yz 






xpar (71-73): Values of TGCs: Set 3, see CERN 96-01, Vol. 1, p. 525 


alpha_Wphi 


=xpar(71) = q;^0 




alpha_Bphi 


=xpar(72) = as^, 




alpha_W 


=xpar(73) = aw 




amel 


=xpar(100) (=0.510998902d-3): beam (electron) mass 


Alflnvl 


=xpar(101) (=137.03599976) 


1/ctQED at the Thomson limit 


gpicob 


=xpar(101) (=389.379292d6) 


GeV~^ to picobarn translation constant 


BR(1:20) 


xpar (131-139) : W branching ratios (BR); the numbering of entries is: 




1 = nd, 2 = cd, 3 = us, 4 = cs, 


5 = ub, 6 = cb, 7 = evg, 8 = iiVn, 9 = Ti^r 


amaf in(20) 


xpar (500 +10*KF +6): Masses of the W decay products; the used entries KF 




(in PDG notation) are: 




1 = d. 


2 = u, 3 = s, 4 = c, 5 = 6, 6 = t. 




11 = e 


, 12 = fe, 13 = /X, 14 = f^, 15 = T, 16 = l^r- 




Note: The masses of r and v-,- 


have to be independently set to the same 




numerical values in the initialization of TAUDLA 


VCKM(1:3,1:3) 


xpar (111-119): CKM matrix elements (PDG 2000) 



Table 3: The list of input parameters of the YFSWW3 generator in the xpar vector (cent.). 
The default values are in brackets. 
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Parameter 


Position and meaning 




xpar (1011-1013) : Dummy parameters in YFSWW3 (specific to KoralW) 


KeyCul 


xpar(1014) (=2): Coulomb correction switch 




=0 


Coulomb correction is OFF 




=1 


"Standard" Coulomb correction is ON 




=2 


"Screened" Coulomb correction is ON 


KeyBra 


xpar (1021) (=2): Sets W branching ratios (BR), used for normalization of 






the matrix element 




=0 


Born values (no mixing) with "naive" QCD (if as = 0: ud,cs = 1/3, 






eu,^v,Tv = 1/9, others = 0) 




=1 


Values of BRs taken from the input 




=2 


With CKM mixing and "naive" QCD, calculated in the IB A from the 






CKM matrix (PDG 2000) 


KeyMas 


xpar (1022): Dummy parameter in YFSWW3 (specific to KoralW) 


KeyZet 


xpar(1023) (=1) : Z width choice 




=0 


Z width in Z propagator: {s/Mz)Tz 




=1 


Z width in Z propagator: MzTz 




=2 


no Z width in Z propagator 


KevSpn 


xpar (1024) (=1): Spin effects in decays 




=0 


Spin effects OFF (for tests only) 




=1 


Spin effects ON (recommended) 


KeyRed 


xpar (1025): Dummy paremeter in YFSWW3 (specific to KoralW) 


KeyWu 


xpar(1026) (=1) : width choice 




=0 


W width in W propagator: {s/M]y)rw 




=1 


W width in W propagator: M\yTw 




=2 


No W width {Tw = 0) in 11^ propagator (dangerous!) 


KevWet 


xpar (1031) (=0): Unweighted/Weighted choice 




=0 


Unweighted events WtMod=l, for apparatus Monte Carlo 




=1 


Weighted events (variable weights) 


KeyRnd 


xpar (1032) (=1): Random number generator switch 




=1 


RANMAR random number generator 




=2 


ECURAN random number generator 




=3 


CARRAN random number generator 


KeySmp 


xpar (1033): Dummy paremeter in YFSWW3 (specific to KoralW) 


KeyMix 


xpar (1041) (=1): Input parameter scheme (IPS) switch 




=0 


LEP2 '95 Workshop scheme (only for Born and ISR) 




=1 


G^-scheme (recommended; the only choice for 0{a) EW corrections) 




=2 


a-scheme ("pure" Born; for tests) 


Key4f 


xpar (1042): Dummy parameter in YFSWW3 (specific to KoralW) 



Table 4: The list of input parameters of the YFSWW3 generator in the xpar vector (cont.). 
The default values are in brackets. 
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Parameter 



Position and meaning 



KeyAcc 



KeyZon 
KeyWon 
KeyDwm 



KeyDwp 
Nout 
Jakl 
Jak2 



Itdkrc 



IfPhot 



IfHadM 
IfHadP 



Umask 

BE params 



m 



W+ 



xpar(1043) (=0): Anomalous TGCs switch 
=0: Anomalous WWV couplings OFF 
>0: Anomalous WWV couplings ON 

=1: The most general (complex number) TGCs in the notation of Hagiwara 

et al., Nucl. Phys. B282 (1987) 253 
=2 Parametrization of CERN 96-01, Vol. 1, p. 525: Sz^x-y^xz^y-^^yz 
=3: Parametrization of CERN 96-01, Vol. 1, p. 525: aw(i>-,oiB(j)-,oiw 
xpar(1044): Dummy parameter in YFSWW3 (specific to KoralW) 
xpar(1045): Dummy parameter in YFSWW3 (specific to KoralW) 
xpar(1055) (=0): Sets decay channel of W~ 
=0: Inclusive; otherwise the exclusive modes for the W are: 
= 12 34567 8 9" 
ud cd us cs ub cb eu iiv tv 
xpar(1056) (=0): Sets the decay channel of ^ similar assignments as for 

w- 

(=-1): Output unit number for the generator (if < then 
Nout=16) 

(=0): Input for TAUOLA, defines decay mode of the r+ 

decay 

(=0): Input for TAUOLA, defines decay mode of the r" 
decay 

Jakl,Jak2 = -1: TAUOLA is switched OFF 

JaJsl , JaJs2 = 0: Requests all decay channels to be simulated 
Jakl, Jak2 > 0: Single specific decay channel, see TAUOLA manual 
xpar (1073) (=1) : Input for TAUOLA, radiative corrections in leptonic r decays 
switch 

Itdkrc = 1: Corrections are ON 

Itdkrc = 0: Corrections are OFF 

xpar (1074) (=1): PHOTOS activation switch 

=0: PHOTOS is OFF 

=1: Radiation in leptonic W decays is ON 

=2: Radiation in both leptonic and quarkonic W decays is ON (for tests) 
xpar (1075) (=1): W~ hadronization activation switch (JETSET) 
xpar (1075) (=1): hadronization activation switch (JETSET) 
IfHadM, if hadP=0: Hadronization is OFF 
IfHadM, if hadP=l: Hadronization is ON 

In the present version IfHadM and IfHadP have to be equal! 
xpar (1101-1302): Dummy parameters in YFSWW3 (specific to KoralW) 
xpar (4061-4070): Dummy parameters in YFSWW3 (specific to KoralW) 



xpar (1057) 
xpar(1071) 
xpar (1072) 



in W- 



Table 5: The list of input parameters of the YFSWW3 generator in the xpar vector (cont.). 
The default values are in brackets. 



40 



Parameter 


Meaning 


Ql(4) 
Q2(4) 

SPhum(4) 
SPhot (100,4) 
Nphot 


Four-momentum of W~ 
Four-momentum of W'^ 

Sum of four-momenta of photons from the WW production stage 
Four-momenta of photons from the WW production stage 
Number of photons from the WW production stage 



Table 6: The list of four-momenta in the COMMON /MOMWWP/ of the YFSWW3 generator. They 
are given in GeV in the CMS of the incoming beams. 



Parameter 


Meaning 


pfl(4) 
pf2(4) 

pf3(4) 
pf4(4) 
Phot (100,4) 
Npho 


Four-momentum of the fermion from the W~ decay 
Four-momentum of the antifermion from the W~ decay 
Four-momentum of the fermion from the W~^ decay 
Four-momentum of the antifermion from the W~^ decay 
Four-momenta of photons from the WW production stage 
Multiphcity of photons from the WW production stage 



Table 7: The list of four-momenta in the COMMON /MOMDEC/ of the YFSWW3 generator. They 
are given in GeV in the CMS of the incoming beams. 



Routine 




Meaning 


YFSWW.GetBeams 


ql(4) 
q2(4) 


Four-momentum of the e~ beam 
Four-momentum of the e"*" beam 


YFSWW_Get4f 


flav(4) 

pi (4) 

p2(4) 
p3(4) 
p4(4) 


Flavours of the final 4/-state in the PDG conven- 
tion 

Four- momentum of the fermion of f lav(l) 
Four-momentum of the fermion of flav(2) 
Four-momentum of the fermion of f lav(3) 
Four-momentum of the fermion of flav(4) 


YFSWW.GetPhotAlI 


NphAll 

PhoAll(100,4) 


Number of photons from the WW production stage 
Four-momenta of these photons 



Table 8: The list of output parameters of the YFSWW3 generator transferred through parameters 
of the getter-routines: YFSWW_Get Beams (ql,q2), YFSWW_Get4f (flav,pl,p2,p3,p4) and 
YFSWW_Get Phot All (NphAll , PhoAll) . 
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Parameter 


Meaning 


WtCrul 
WtCru2 

WtMod 


Crude weight for ISR 
Crude weight for the rest 

They are necessary to build the total weight out of WtSet 
The "best" weight corresponding to Eq. (^), e.g. for KeyCor=5: 
WtMod = Wtcrul*Wtcru2*(WtSet(41)-WtSet(2)+WtSet(4)) 


WtSet(l-19) 

WtSet(l) 

WtSet(2) 

WtSet(3) 

WtSet(4) 

WtSet(lO) 

WtSet(ll) 

WtSet(12) 

WtSet(13) 

WtSet(14) 

WtSet(15) 

WtSet(16) 

WtSet(17) 

WtSetdS) 

WtSet(19) 


Various ISR-type contributions to the matrix element {(3 func- 
tions) - to get the total weight they must be multiplied by WtCrul 
Zeroth-order (/3o) 
First-order ISR (/3o + h) 
becond-order ibR (po + Pi + P2J 
Third-order ISR (^0 + + 
0[a^) contribution from /3o 
0{a^) ISR contribution from /3o 
0{a^) ISR contribution from 
O(a^) ISR contribution from 
O(a^) ISR contribution from /3i 
O(a^) ISR contribution from ^2 
0\a ) iSR contribution from po 
Oic?^ ISR contribution from ^\ 
O(a^) ISR contribution from ^2 
0{g?^ ISR contribution from /^a 


WtSet (41-52) 

WtSet(41) 
WtSet(51) 
WtSet(42) 
WtSet(43) 
WtSet(52) 


Contributions to the matrix element from 0{ol) EW corrections 
m the vv W production stage ~ to get the total weight they must 
be multiphed by WtCrul *WtCru2 
0{a^) EW corrections (/3o + 

Pretabulated ("fast") 0{a^) EW corrections (^0 + A) 
0{a^) contribution from Po 
0{a^) contribution from /3i 
Pretabulated 0{a^) contribution from (]q 


WtSet (65-68) 
WtSet(65) 
WtSet (66) 
WtSet (67) 
WtSet (68) 


Fixed-order (no exponentiation) contributions (for tests) 
0{a^), no exponentiation 
0{a^), no exponentiation 

0{a^) virtual -|- soft real photon only, no exponentiation 
0{a^) real hard photon only, no exponentiation 



Table 9: The list of output weights in the COMMON /WGTALL/ of the YFSWW3 generator. 
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Routine 


Parameter 


Meaning 


YFSWW_GetWtMain 


WtMain 


Principal best MC event weight 


YFSWW_GetWtISR 


WtlSR 


Best ISR-type weight 


YFSWW_GetWtAll 


WtMain 
WtCrud 

WtSetAll(lOO) 


Principal best MC event weight 

Crude MC event weight: WtCrud=WtCrul*WtCru2 

Model weights 

(the same as in WtSet(lOO) in COMMON /WGTALL/) 



Table 10: The list of output parameters of the YFSWW3 generator transferred through pa- 
rameters of the getter-routines: YFSWW_GetWtMain (WtMain) , YFSWW_GetWtISR(WtISR) and 
YFSWW_GetWt All (WtMain , WtCrud , WtSet All) . 



Routine 


Parameter 


Meaning 


YFSWW_GetXSecMC 


XSecMC 
XErrMC 


Principal total MC cross section [pb] 
Its absolute error [pb] 


YFSWW_GetXSecNR 


XSecNR 
XErrNR 


Normalization (primary) cross section, see Sect. ^^: 
For KeyWgt=0, the principal cross section XSecMC [pb], 
For KeyWgt=l, the crude cross section XCrude [pb] 
Its absolute error [pb] 


YFSWW_GetNevMC 


NevMC 


Total number of generated Monte Carlo events 



Table 11: The list of output parameters of the YFSWW3 generator transferred 
through parameters of the getter-routines: YFSWW_GetXSecMC (XSecMC, XErrMC) , 
YFSWW_GetXSecNR (XSecNR, XErrNR) and YFSWW_GetNevMC (NevMC) . 
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B Output of the Demo Program 



* YFSWW3 version 1.16. January 2001 * 

* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 



YFSWW3 version 1.16, January 2001 
S. Jadach, W. Placzek, M. Skrzypek, 
B.F.L. Ward, Z. Was 

INPUT PARAMATERS: 

200.00000000 CMS energy total 

5 Rad. Corr. switch 

2101 QED super-switch 

1 Init . state Rad . 

Final state Rad. 

1 Next-To-Leading 

2 Coulomb corr. 

101112 Physics super-switc 

FS mass reduction 

1 Spin in W decays 
1 Z propag. 

1 Mass kinematics . 

2 Branching Rat . 
1 W propag. 

10 Technical super-swi 

LPA_a or LPA_b mode 

Presampler type 

1 Rand. Numb. Gen. 

Weighting switch 

10001 Miscelaneous 

1 EW Input Par. Schem 
4 fermion matr el 

Anomalous couplings 

1 WW type final state 
ZZ type final state 

W- decay mode 

W+ decay mode 



1.16639000 
132.50494581 

91.18820000 
2.49520000 

80.41900000 
2.09957845 
0.00000100 
0.99000000 
4.00000000 
4.00000000 
0.11850000 



G_mu * Ids 
Inverse alpha_w 

Z mass [GeV] 
Z width [GeV] 
W mass [GeV] 
W width [GeV] 
Dummy infrared cut 
v_max ( =1 ) 
Max wt for rejectn. 
Max wt for CC03 rej 
alpha_s : QCD coupl . 
Color Re-Con. Prob. 



*********************************************** 

*******i****f **f *****************|**1**|*|*****^ 
******* *********** * * *********** ******* ******** * 



0.97493000 
0.22250000 
0.00350000 
-0.22246000 



CKM matrix elements : 



V_ud 
V_us 
V_ub 
V_cd 



CMSEne 

KeyCor 

KeyRad 
KeylSR 
KeyFSR 
KeyNLL 
KeyCul 

KeyPhy 
KeyRed 
KeySpn 
KeyZet 
KeyMas 
KeyBra 
KeyWu 

KeyTek 
KeyLPA 
KeySmp 
KeyRnd 
KeyWgt 

KeyMis 

KeyMix 

Key4f 

KeyAcc 

KeyWon 

KeyZon 

KeyDWm 
KeyDWp 

Gmu 

alfWin 

aMaZ 

GammZ 

aMaW 

GammW 

VVmin 

VVmax 

WtMax 

WtMax 

alphas 

PReco 

sinW2 



VCKM(1,1) 
VCKM(1,2) 
VCKM(1,3) 
VCKM(2,1) 



I.O 

IRl 

IQl 
IQ2 
IQ3 
IQ4 
IQ5 

IPl 
IP2 
IP3 
IP4 
IP5 
IP6 
IP7 

ITl 

IT2 
ITS 
IT4 
IT5 

IMl 
IM2 
IM3 
IM4 
IM5 
IM6 

IDl 
ID2 

I.l 
1.2 
1.3 



110 
111 
112 

113 



IVl 
IV2 
IV3 
IV4 
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* 





97412000 




V 


cs 


VCKM(2,2) 


IV5 


* 


* 





04000000 




V 


cb 


VCKM(2,3) 


IV6 


* 


* 





00549000 




V 


.td 


VCKM(3,1) 


IV7 


* 


* 


-0 


03978000 




V 


ts 


VCKM(3,2) 


IV8 


* 


* 





99920000 




V. 


tb 


VCKM(3,3) 


IV9 


* 


* 




Unitarity check of 


the CKM 


matrix: 






* 


* 




1.00001 


0.00000 


0.00000 








* 


* 




VV+ = 0.00000 


1.00000 


0.00000 








* 


* 
* 




0.00000 


0.00000 


1.00001 








* 
* 



1 

4 
174 






0.32071536 
0.01669847 
0.01670448 
0.32018266 
0.00000413 
0.00053987 
0.10838559 
0.10838559 
0.10838559 



0.00600000 
0.00300000 
0.12250000 
25000000 
20000000 
30000000 
00051100 
00000000 
10565836 
0.00000000 
1.77703000 
0.00000000 



Branching ratios of W decays: 



115.00000000 





1 
1 
1 
1 



fermion masses: 



* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

* _ 



Higgs mass : 



ud 


BR(1) 


IBl 


* 


cd 


BR(2) 


IB2 


* 


us 


BR(3) 


IB3 


* 


cs 


BR(4) 


IB4 


* 


ub 


BR(5) 


IB5 


* 


cb 


BR(6) 


IB6 


* 


e 


BR(7) 


IB7 


* 


mu 


BR(8) 


IB8 


* 


tau 


BR(9) 


IB9 


* 
* 


d 


amaf in ( 1 ) 


IMl 


* 


u 


amaf in(2) 


IM2 


* 


s 


amaf in (3) 


IMS 


* 


c 


amaf in (4) 


IM4 


* 


b 


amaf in (5) 


IM5 


* 


t 


amaf in(6) 


IMS 


* 


e 


amaf in ( 1 1 ) 


IM7 


* 


ve 


amaf in(12) 


IMS 


* 


mu 


amaf in(13) 


IM9 


* 


vmu 


amaf in (14) 


IMIO 


* 


tau 


amaf in(15) 


IMll 


* 


vtau 


amaf in (16) 


IM12 


* 
* 


Higgs 


aMHig 


IMH 


* 



DECAY LIBRARIES * 

TAUOLA for W+ Jakl ILl * 

TAUOLA for W- Jak2 IL2 * 

TAUOLA Ord (alpha) Itdkrc ILS * 

PHOTOS IfPhot IL4 * 

JETSET for W- IfHadM IL5 * 

JETSET for W+ IfHadP IL6 * 



***###***###*######^^ 

### ### ######## ######## ### 



* 
* 
* 
* 
* 
* 
* 



* ### ### ######## ######## ### ### ### 

* ### ### ### #### ### ### ### 

* ### ###### #### ### ### ### 

* ### ###### #### ### # ### ### 

* ### ### ######## #### #### #### 

* ### ### ###### ### ### ### 



### 
### 
### 
### 
# ### 
#### 
### 



##### 
####### 

### 
### 
### 
####### 

##### 



* Initialize KarLud start * 

* 19.33542764 xs.crude VESKO XCVESK * 

* 19.19161478 xs_crude GAUSS XCGAUS * 

* 0.00749353 XCVESK/XCGAUS-1 * 



* Initialize KarFin start * 



45 



* *****TAUOLA LIBRARY: VERSION 2.6 ****** * 

* ***********Augus"t 1995*************** * 

* **AUTHDRS: S.JADACH, Z . WAS************* * 

* **R. DECKER, M. JEZABEK, J . H . KUEHN***** * 

* **AVAILABLE FROM: WASM AT CERNVM ****** * 

* ***** PUBLISHED IN COMP . PHYS . COMM.*** * 

* *******CERN-TH-5856 SEPTEMBER 1990***** * 

* *******CERN-TH-6195 SEPTEMBER 1991***** * 

* *******CERN-TH-6793 NOVEMBER 1992***** * 

* **5 or more pi dec: precision limited * 

* ******DEXAY ROUTINE: INITIALIZATION**** * 

* JAKl = DECAY MODE FERMIONl (TAU+) * 

* JAK2 = DECAY MODE FERMI0N2 (TAU-) * 



YFSWW demo 

NevTot = 1000 events to be generated . . . 



skipped .... 

Event listing (standard) 



I 


particle/jet 


K(I,1) 


K(I,2) 


K(I,3) 




K(I,4) 




K(I,5) 


Pd.i) 


P(I,2) 


P(I,3) 


P(I,4) 


P(1.5) 


1 


!e-l 




21 


11 









3 






4 





00000 





00000 


100 


00000 


100 


00000 





00061 


2 


>e-f ! 




21 


-11 









3 






4 





00000 





00000 


-100 


00000 


100 


00000 





00061 


3 


(M-) 




11 


-24 


1 






B 






6 


-68 


16061 





61640 


33 


01304 


96 


19679 


67 


76468 


4 


(M+) 




11 


24 


1 






7 






8 


58 


15061 


-0 


61640 


-33 


01304 


104 


80321 


80 


69684 


6 


(d) 




14 


1 


3 


3 


6 


10 








10 


-62 


20078 


-14 


63227 


6 


34668 


64 


21308 





00600 


6 


(u-) 




14 


-2 


3 








11 


3 


5 


11 


4 


05017 


15 


24767 


26 


66636 


30 


98371 





00300 


7 


(u) 




14 


2 


4 


3 


8 


29 








29 


18 


03726 


6 


62502 


26 


47908 


32 


69649 





00300 


8 


(d-) 




14 


-1 


4 








30 


3 


7 


30 


40 


11336 


-7 


14042 


-69 


49212 


72 


10672 





00600 


9 


CCMshower) 




11 


94 


6 






10 






11 


-58 


15062 





61640 


33 


01304 


96 


19679 


67 


76468 


10 


(d) 




14 


1 


9 


3 


5 


13 





5 


12 


-60 


91610 


-14 


19066 


6 


47966 


63 


22975 


6 


62341 


11 


(n-) 




14 


-2 


9 





6 


14 


3 


6 


15 


2 


76548 


14 


80606 


26 


53338 


31 


96705 


9 


53966 


12 


(d) 




14 


1 


10 


3 


13 


17 





10 


16 


-46 


33492 


-8 


95829 


6 


22708 


47 


67010 


2 


54677 


13 


(g) 




14 


21 


10 


3 


10 


18 


3 


12 


19 


-14.68117 


-6 


23237 





25269 


16 


65966 


1 


43197 


14 


(u-) 






-2 


11 









3 






2 


79617 
















90286 


15 


(g) 




13 


21 


11 


2 


14 





2 


11 





-0 


02970 


-0 


14691 





54657 





56675 





00000 


16 


Cd) 




13 


1 


12 


2 


17 








12 





-12 


64723 


-1 


30492 


1 


82179 


12 


84423 





00990 


17 


Cg) 




13 


21 


12 


2 


12 





2 


16 





-33 


68769 


-7 


65337 


4 


40528 


34 


82587 





00000 


18 


Cg) 




13 


21 


13 


2 


13 





2 


19 





-9 


68361 


-4 


03725 


-0 


24096 


10 


49428 





00000 


19 


Cg) 




13 


21 


13 


2 


18 





2 


13 





-4 


89756 


-1 


19512 





49355 


5 


06537 





00000 


20 


(U-) 




14 


-2 


14 





14 


22 


3 


21 


23 


4 


17251 


13 


40366 


18 


62694 


23 


57434 


3 


42331 


21 


Cg) 




14 


21 


14 


3 


20 


24 


3 


14 


26 


-1 


37734 


1 


64930 


7 


36988 


7 


82596 


1 


66749 


22 


(u-) 




14 


-2 


20 





20 


26 


3 


23 


27 


4 


29870 


12 


69084 


16 


24373 


21 


06317 


1 


63932 


23 


(g) 




13 


21 


20 


2 


22 





2 


20 





-0 


12619 





81281 


2 


38321 


2 


52117 





00000 


24 


Cg) 




13 


21 


21 


2 


21 





2 


25 





-1 


25460 





15914 


3 


44828 


3 


67287 





00000 


25 


(g) 




13 


21 


21 


2 


24 





2 


21 





-0 


12274 


1 


39016 


3 


91159 


4 


15309 





00000 


26 


(u-) 




13 


-2 


22 





22 





2 


27 





2 


89930 


7 


60867 


8 


72410 


11 


86990 





00660 


27 


(g> 




13 


21 


22 


2 


26 





2 


22 





1 


39939 


6 


08227 


7 


61962 


9 


18327 





00000 


28 


CCMshower) 




11 


94 


7 






29 






30 


58 


15061 


-0 


61540 


-33 


01304 


104 


80321 


80 


69683 


29 


(u) 




14 


2 


28 


3 


7 


32 





7 


31 


23 


86922 


5 


09816 


15 


70427 


43 


16896 


31 


95629 


30 


(d-) 




14 


-1 


28 





8 


33 


3 


8 


34 


34 


28139 


-5 


71356 


-48 


71731 


61 


63424 


14 


74939 


31 


(u) 




13 


2 


29 


2 


32 








29 





2 


06978 


-3 


06418 


-6 


61677 


6 


71930 





00560 


32 


(g) 




14 


21 


29 


3 


29 


35 


3 


31 


36 


21 


79944 


8 


16234 


21 


32004 


36 


44967 


18 


23068 


33 


(d-) 




14 


-1 


30 





30 


37 


3 


34 


38 


9 


91990 


3 


71419 


-14 


20968 


17 


89177 


2 


46019 


34 


Cg) 




14 


21 


30 


3 


33 


39 


3 


30 


40 


24.36149 


-9 


42776 


-34 


60773 


43 


74247 


6 


34476 


35 


(g) 




14 


21 


32 


3 


32 


41 


3 


36 


42 


24 


19235 


8 


81966 


20 


09323 


33 


02196 


4 


86387 


36 


(g) 




14 


21 


32 


3 


35 


43 


3 


32 


44 


-2 


39291 


-0 


66722 


1 


22681 


3 


42771 


2 


01816 


37 


(d-) 




13 


-1 


33 





33 





2 


38 





7 


66346 


2 


82714 


-12 


47236 


14 


90909 





00990 


38 


(g) 




13 


21 


33 


2 


37 





2 


33 





2 


25644 





88705 


-1 


73723 


2 


98268 





00000 


39 


(g) 




14 


21 


34 


3 


34 


46 


3 


40 


46 


11 


13067 


-1 


66836 


-12 


38914 


16 


78682 


1 


29097 


40 


(g) 




14 


21 


34 


3 


39 


47 


3 


34 


48 


13 


23092 


-7 


76938 


-22 


11869 


26 


96666 


1 


39821 


41 


(g) 




13 


21 


35 


2 


35 





2 


42 





6 


51083 


2 


61590 


3 


06918 


7 


65857 





00000 


42 


(g) 




14 


21 


35 


3 


41 


49 


3 


35 


50 


17 


68152 


6 


20366 


17 


02405 


25 


36338 


1 


53671 


43 


(g) 




13 


21 


36 


2 


36 





2 


44 





-2 


15456 


-0 


80571 





13785 


2 


30441 





00000 


44 


(g) 




13 


21 


36 


2 


43 





2 


36 





-0 


23836 





13849 


1 


08896 


1 


12330 





00000 


46 


(g) 




13 


21 


39 


2 


39 





2 


46 





9 


64866 


-1 


43393 


-9 


90045 


13 


82938 





00000 


46 


(g) 




13 


21 


39 


2 


45 





2 


39 





1 


68193 


-0 


22443 


-2 


48869 


2 


95744 





00000 


47 


(g) 




13 


21 


40 


2 


40 





2 


48 





1 


41318 


-1 


20728 


-3 


26543 


3 


75734 





00000 


48 


(g) 




13 


21 


40 


2 


47 





2 


40 





11 


81774 


-6 


56210 


-18 


85316 


23 


19831 





00000 


49 


(g) 




13 


21 


42 


2 


42 





2 


50 





14 


75801 


5 


72986 


14 


60889 


21 


54181 





00000 


60 


(g) 




13 


21 


42 


2 


49 





2 


42 





2 


92351 





47380 


2 


41516 


3 


82157 





00000 


61 


(d) 


A 


12 


1 


16 






73 






73 


-12 


64723 


-1 


30492 


1 


82179 


12 


84423 





00990 


62 


(g) 


I 


12 


21 


17 






73 






73 


-33 


68769 


-7 


66337 


4 


40628 


34 


82687 





00000 


63 


(g) 


I 


12 


21 


19 






73 






73 


-4 


89756 


-1 


19612 





49366 


6 


06537 





00000 


54 


(g) 


I 


12 


21 


18 






73 






73 


-9 


68361 


-4 


03725 


-0 


24096 


10 


49428 





00000 


55 


(g) 


I 


12 


21 


15 






73 






73 


-0 


02970 


-0 


14691 





54657 





56675 





00000 


56 


(g) 


I 


12 


21 


25 






73 






73 


-0 


12274 


1 


39016 


3 


91159 


4 


15309 





00000 


67 


(g) 


I 


12 


21 


24 






73 






73 


-1 


26460 





16914 


3 


44828 


3 


67287 





00000 


68 


(g) 


I 


12 


21 


23 






73 






73 


-0 


12619 





81281 


2 


38321 


2 


62117 





00000 



46 



59 


(g) 


I 12 


21 


27 


73 


73 


60 


(u-) 


V 11 


-2 


26 


73 


73 


61 


(u) 


A 12 


2 


31 


91 


91 


62 


(g) 


I 12 


21 


44 


91 


91 


63 


(g) 


I 12 


21 


43 


91 


91 


64 


(g) 


I 12 


21 


50 


91 


91 


66 


(g) 


I 12 


21 


49 


91 


91 


66 


(g) 


I 12 


21 


41 


91 


91 


67 


(g) 


I 12 


21 


48 


91 


91 


68 


(g) 


I 12 


21 


47 


91 


91 


69 


(g) 


I 12 


21 


46 


91 


91 


70 


(g) 


I 12 


21 


45 


91 


91 


71 


Cg) 


I 12 


21 


38 


91 


91 


72 


(d-) 


V 11 


-1 


37 


91 


91 


73 


(string) 


11 


92 


61 


74 


90 


74 


(eta) 


11 


221 


73 


123 


125 


75 


nO 


1 


2112 


73 








76 


P"- 


1 


-2212 


73 








77 


(rhoO) 


11 


113 


73 


126 


127 


78 


K+ 


1 


321 


73 








79 


(LafflbdaO) 


11 


3122 


73 


128 


129 


80 


n-O 


1 


-2112 


73 








81 


(eta') 


11 


331 


73 


130 


131 


82 


(piO) 


11 


111 


73 


132 


133 


83 


(omega) 


11 


223 


73 


134 


136 


84 


(piO) 


11 


111 


73 


137 


138 


85 


(rhoO) 


11 


113 


73 


139 


140 


86 


pi- 


1 


-211 


73 








87 


(rho+) 


11 


213 


73 


141 


142 


88 


pi- 


1 


-211 


73 








89 


K+ 


1 


321 


73 








90 


(K*-) 


11 


-323 


73 


143 


144 


91 


(string) 


11 


92 


61 


92 


122 


92 


(rho+) 


11 


213 


91 


146 


146 


93 


pi- 


1 


-211 


91 








94 


pit 


1 


211 


91 








95 


(rhoO) 


11 


113 


91 


147 


148 


96 


(rhoO) 


11 


113 


91 


149 


150 


97 


(rho-) 


11 


-213 


91 


151 


152 


98 


(Delta++) 


11 


2224 


91 


153 


154 


99 


K- 


1 


-321 


91 








100 


(Sigma~-) 


11 


-3222 


91 


155 


156 


101 


P+ 


1 


2212 


91 








102 


(rho+) 


11 


213 


91 


157 


158 


103 


P'- 


1 


-2212 


91 








104 


(eta) 


11 


221 


91 


169 


161 


105 


(rho-) 


11 


-213 


91 


162 


163 


106 


(rho+) 


11 


213 


91 


164 


165 


107 


(eta) 


11 


221 


91 


166 


168 


108 


(piO) 


11 


111 


91 


169 


170 


109 


(rho-) 


11 


-213 


91 


171 


172 


110 


(rho+) 


11 


213 


91 


173 


174 


111 


(DeltaO) 


11 


2114 


91 


175 


176 


112 


Ceta) 


11 


221 


91 


177 


179 


113 


(Delta--) 


11 


-2214 


91 


180 


181 


114 


(eta) 


11 


221 


91 


182 


183 


115 


pit 


1 


211 


91 








116 


(piO) 


11 


111 


91 


184 


185 


117 


(rho-) 


11 


-213 


91 


186 


187 


118 


(Sigma*+) 


11 


3224 


91 


188 


189 


119 


(Xi»-0) 


11 


-3324 


91 


190 


191 


120 


K- 


1 


-321 


91 








121 


(rhoO) 


11 


113 


91 


192 


193 


122 


pi+ 


1 


211 


91 








123 


(piO) 


11 


111 


74 


194 


196 


124 


(piO) 


11 


111 


74 


196 


197 


125 


(piO) 


11 


111 


74 


198 


199 


126 


pi+ 


1 


211 


77 








127 


pi- 


1 


-211 


77 








128 


nO 


1 


2112 


79 








129 


(piO) 


11 


111 


79 


200 


201 


130 


ganuna 


1 


22 


81 








131 


(rhoO) 


11 


113 


81 


202 


203 


132 


gamma 


1 


22 


82 








133 


gamma 


1 


22 


82 








134 


pi- 


1 


-211 


83 








13B 


pit 


1 


211 


83 








136 


(piO) 


11 


111 


83 


204 


205 


137 


gamma 


1 


22 


84 








138 


gamma 


1 


22 


84 








139 


pi- 


1 


-211 


85 








140 


pit 


1 


211 


86 








141 


pi+ 


1 


211 


87 








142 


(piO) 


11 


111 


87 


206 


207 


143 


(K-O) 


11 


-311 


90 


208 


208 


144 


pi- 


1 


-211 


90 








146 


pit 


1 


211 


92 








146 


(piO) 


11 


111 


92 


209 


210 



1 


.39939 


5 


.08227 


7 


.51962 


9 


. 18327 





.00000 


2 


.89930 


7 


.60867 


8 


.72410 


11 


.86990 





.00660 


2 


.06978 


-3 


.06418 


-6 


.61677 


6 


.71930 





.00660 


-0 


. 23835 





. 13849 


1 


.08896 


1 


. 12330 





. 00000 


-2 


. 15456 


-0 


.80571 





. 13785 


2 


.30441 





. 00000 


2 


.92351 





.47380 


2 


.41516 


3 


.82157 





. 00000 


14 


.75801 


5 


.72986 


14 


. 60889 


21 


.54181 





. 00000 


6 


.51083 


2 


.61590 


3 


.06918 


7 


.65857 





. 00000 


11. 


.81774 


-6 


.56210 


-18. 


.85316 


23 


.19831 





.00000 


1 


.41318 


-1 


.20728 


-3. 


.26643 


3 


.75734 





.00000 


1 


.58193 


-0 


.22443 


-2 


.48869 


2 


.95744 





.00000 


9 


.54865 


-1 


.43393 


-9 


.90045 


13 


. 82938 





. 00000 


2 


.25644 





.88705 


-1 


.73723 


2 


.98268 





. 00000 


7 


.66346 


2 


.82714 


-12. 


.47235 


14 


.90909 





.00990 


-68, 


. 16062 





.61639 


33 


.01306 


96 


.19680 


67 


.76468 


-10. 


. 24036 


-1 


.30218 


1 


.37509 


10 


.42839 





.64760 


-18 


. 93694 


-4 


. 18797 


2 


.49748 


19 


.57721 





. 93960 


-14 


.20916 


-3 


.08438 


1 


.25868 


14 


. 62458 





. 93830 


-6. 


.22881 


-2 


.31614 


1 


.02168 


6 


.76712 





.76646 


-3, 


.92216 


-0 


.97100 





.06270 


4 


.07107 





.49360 


-3, 


.01880 


-0 


.44703 





.30808 


3 


.26381 


1 


.11660 


-3. 


.76048 


-1 


.33794 





.39686 


4 


.11053 





.93960 


-0. 


. 12766 





.40474 





.86328 


1 


.35133 





.96807 


-0 


.50601 


-0 


.14221 





.41852 





.68531 





. 13600 


-0 


. 63008 





.00890 


3 


.96938 


4 


. 09435 





.78146 





.12911 





.68122 





.54176 





. 89020 





. 13500 





.01639 


1 


.31382 


3 


.53153 


3 


.84545 





.76771 


-0, 


.26161 


-0 


.09021 





.26972 





.39798 





.13960 





.28192 


1 


.27622 


2 


.20383 


2 


.67551 





.77191 


-0 


.03182 


1 


.50185 


1 


.71206 


2 


.28193 





. 13960 


1 


.10097 


1 


.59131 


2 


.64052 


3 


.31065 





. 49360 


2 


. 17486 


7 


.71740 


9 


.96190 


12 


.82138 





.92739 


68, 


. 16061 


-0 


.61640 


-33. 


.01306 


104. 


.80321 


80 


.69683 


1 


.65024 


-2 


.13893 


-3 


. 14091 


4 


.17312 





.76686 


-0 


.17159 





.04685 





.08614 





.24196 





. 13960 





.24570 


-0 


.03059 


-0 


.39585 





.48733 
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* window A * 

* 4537 total no of events NEVTOT AO * 

* WT<0 events NEVNEG Al * 

* 19.33542764 xs.cru VESKO [pb] XCVESK A2 * 

* 19.17361574 +- 0.01754039 xs_est VESKO [pb] XSVE A3 * 

* 19.19161478 +- 0.00001919 xs_est Gauss [pb] XSGS A4 * 

* -0.00093786 +- 0.00091582 XCVE/XCGS-1 A5 * 

* KarFin output - window A * 

* Weight Statistics * 

* 0.99643450 +- 0.00368336 general weight WT Al * 

* 4263 generated events NEVGEN A2 * 

* 0.08139808 aver. ph. multi. AVMULT A3 * 

* Marked photons MARTOT A4 * 

* 0.99108609 +- 0.00145252 Kinematics, smin WTl A5 * 

* 0.99999743 +- 0.00000168 Jacobian WT2 A6 * 

* 1.00534843 +- 0.00335486 Photon ang. dist . WT3 A7 * 

* KarFin output - window B * 

* on mass weights * 

* 0.96568167 +- 0.00245351 removal wgt WTREM WT6 Bl * 

* 4263 no . of raw events B2 * 

* WT6=0 events B3 * 

* 0.08139808 raw ph. multipl. B4 * 

* 0.99951442 +- 0.00245437 control wgt WCTRL WT5 B5 * 

* 0. lOOOOOOOE-05 epsilon B6 * 

* 0. lOOOOOOOE-08 delta B7 * 



* ******** DADMEL FINAL REPORT ******** * 

* 48 NEVRAW = NO . OF EL DECAYS TOTAL * 

* 18 NEVACC = NO. OF EL DECS. ACCEPTED * 

* NEVOVR = NO. OF OVERWEIGHTED EVENTS * 

* 0.48927E-12 PARTIAL WTDTH ( ELECTRON) IN GEV UNITS * 

* 1.207993507 IN UNITS GFERMI**2*MASS**5/192/PI**3 * 

* 0.128391966 RELATIVE ERROR OF PARTIAL WIDTH * 

* COMPLETE QED CORRECTIONS INCLUDED * 

**********************************************^ 

skipped .... 

*******************************************^ 
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============ demo ============ 

16.96445667 +- 0.24052939 MC Best, XPAR, YFSWW3 

18.96537427 +- 0.00085125 SemiAnal Born, KorWan 

17.42293229 +- 0.00255375 SemiAnal 0(alf3)exp.LL, KorWan 

========== End demo ========== 



51 



References 

[1] Physics at LEP2, edited by G. Altarelli, T. Sjostrand, and F. Zwirner (CERN 96-01, 
Geneva, 1996), 2 vols. 

[2] Reports of the Working Groups on Precision Calculations for LEP2 Physics, edited 
by S. Jadach, G. Passarino, and R. Pittau (CERN 2000-009, Geneva, 2000). 

[3] S. Jadach, W. Placzek, M. Skrzypek, and B. F. L. Ward, Phys. Rev. D54, 5434 
(1996). 

[4] S. Jadach et al, Phys. Lett. B417, 326 (1998). 
[5] S. Jadach et al, Phys. Rev. D61, 113010 (2000). 

[6] S. Jadach et al., Precision Predictions for (Un) Stable W~^W~ Production At and 
Beyond LEP2 Energies, preprint CERN-TH/2000-337; [hep-ph/OOOTOll ; to be sub- 
mitted to Phys. Lett. B. 

[7] J. Fleischer, F. Jegerlehner, and M. Zralek, Z. Phys. C42, 409 (1989). 

[8] K. Kolodziej and M. Zralek, Phys. Rev. D43, 3619 (1991). 

[9] J. Fleischer, K. Kolodziej, and F. Jegerlehner, Phys. Rev. D47, 830 (1993). 

[10] J. Fleischer, K. Kolodziej, and F. Jegerlehner, Phys. Rev. D49, 2174 (1994). 

[11] J. Fleischer, F. Jegerlehner, K. Kolodziej, and G. J. van Oldenborgh, Comput. Phys. 
Commun. 85, 29 (1995). 

[12] A. Denner, S. Dittmaier, M. Roth, and D. Wackeroth, Phys. Lett. B475, 127 (2000). 

[13] A. Denner, S. Dittmaier, M. Roth, and D. Wackeroth, Nucl. Phys. B587, 67 (2000). 

[14] W. Beenakker, F. A. Berends, and A. Chapovsky, Nucl. Phys. B548, 3 (1999). 

[15] W. Beenakker et al, WW Cross-Sections and Distributions, in Ref. [Q, Vol. 1, p. 79. 

[16] A. Vicini, Acta Phys. Polon. B29, 2847 (1998). 

[17] D. Bardin et al., Event Generators for WW Physics, in Ref. [|l|. Vol. 2, p. 3. 

[18] M. Griinewald et al., Four-Fermion Production in Electron- Positron Collisions, in 
Ref. 0],p. 1. 

[19] R. G. Stuart, Nucl. Phys. B498, 28 (1997), and references therein. 

[20] W. Placzek et al, W-Paii Production with YFSWW/KoralW, preprint CERN- 
TH/2000-357; |hep-ph/0012094| . 



52 



[21] M. Skrzypek, S. Jadach, W. Placzek, and Z. W§,s, Comput. Phys. Commun. 94, 216 
(1996). 

[22] M. Skrzypek et al, Phys. Lett. B372, 289 (1996). 

[23] S. Jadach et al, Comput. Phys. Commun. 119, 272 (1999). 

[24] M. Skrzypek and Z. W^s, Comput. Phys. Commun. 125, 8 (2000). 

[25] S. Jadach et al., The Monte Carlo Program KoralW version 1.51 and The Con- 
curent Monte Carlo KoralW*YFSWW3 with All Background Graphs and First Order 
Corrections to W^-pair Production, preprint CERN-TH/2001-040, UTHEP-01-0102, 
February 2001, (to be submitted to Comput. Phys. Commun.). 

[26] S. Jadach et al, in preparation. 

[27] E. Barbcrio, B. van Eijk, and Z. W§s, Comput. Phys. Commun. 66, 115 (1991), ibid. 
79 291 (1994). 

[28] V. S. Fadin, V. A. Khoze, A. Martin, and W. Stirling, Phys. Lett. B363, 112 (1995). 

[29] A. P. Chapovsky and V. A. Khoze, Eur. Phys. J. C9, 449 (1999). 

[30] W. Bccnakkcr, A. P. Chapovsky, and F. A. Berends, Nucl. Phys. B508, 17 (1997). 

[31] W. Beenakker, A. P. Chapovsky, and F. A. Berends, Phys. Lett. B411, 203 (1997). 

[32] A. Denner, S. Dittmaier, and M. Roth, Nucl. Phys. B519, 39 (1998). 

[33] A. Denner, S. Dittmaier, and M. Roth, Phys. Lett. B429, 145 (1998). 

[34] S. Jadach and B. F. L. Ward, Comput. Phys. Commun. 56, 351 (1990). 

[35] S. Jadach and B. F. L. Ward, Phys. Lett. B274, 470 (1992). 

[36] S. Jadach, B. F. L. Ward, and Z. W^s, Comput. Phys. Commun. 79, 503 (1994). 

[37] B. F. L. Ward, Phys. Rev. D36, 939 (1987), ibid. 42 (1990) 3249. 

[38] K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B274, 1 (1986). 

[39] F. A. Berends et al, Nucl. Phys. B206, 61 (1982). 

[40] T. Kinoshita, J. Math. Phys. 3, 650 (1962). 

[41] T. D. Lee and M. Nauenberg, Phys. Rev. 133, 1549 (1964). 

[42] D. R. Yennie, S. Frautschi, and H. Suura, Ann. Phys. (NY) 13, 379 (1961). 

[43] S. Jadach, E. Richter-Was, B. F. L. Ward, and Z. Was, Comput. Phys. Commun. 
70, 305 (1992). 

53 



[44] S. Jadach, M. Skrzypek, and B. Ward, Phys. Lett. B257, 173 (1991). 

[45] S. Jadach et ai, Comput. Phys. Commun. 102, 229 (1997). 

[46] S. Jadach, W. Placzek, and B. F. L. Ward, Phys. Lett. B390, 298 (1997). 

[47] S. Jadach, B. F. L. Ward, and Z. Wqs, Coheerent Exclusive Exponentiation for 
precision Monte Carlo Calculations, preprint CERN-TH/2000-087; UTHEP-99-0901; 
liep-ph/0006359| ; Phys. Rev. D., in press. 



[48] S. Jadach, B. F. L. Ward, and Z. W§s, Comput. Phys. Commun. 130, 260 (2000). 
[49] S. Jadach, Acta Phys. Polon. B16, 1007 (1985). 



[50] S. Jadach, Practical guide to Monte Carlo, 1999, eprint [physics / 9906056] , also avail- 



able from [http : / /home . cern . ch/ ^ ] adach| (unpublished) . 



[51] G. J. van Oldenborgh and J. A. M. Vermaseren, Z. Phys. C46, 425 (1990). 

[52] R. Decker, S. Jadach, J. H. Kiihn, and Z. W§s, Comput. Phys. Commun. 76, 361 
(1993). 

[53] T. Sjostrand and M. Bengtsson, Comput. Phys. Commun. 43, 367 (1987). 

[54] P. Golonka, E. Richter-W^s, and Z. W§s, The TAUOLA-PHOTOS-F environment 
for versioning the TAUOLA and PHOTOS packages, |hep-ph/0009"302 . 



[55] T. Pierzchala, E. Richter-Wcis, Z. Was, and M. Worek, Spin effects in r-lepton pair 
production at LHC, preprint TPJU 2/2001, January 2001; |hep-ph/OlOim| . 

[56] S. Jadach, The FORTRAN code GlibK, 1995 (unpublished). 



54 



